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ABSTRACT 

We  have  evaluated  the  energetics  of  38  solar  eruptive  events  observed  by  a  variety  of  spacecraft  instruments  between 
2002  February  and  2006  December,  as  accurately  as  the  observations  allow.  The  measured  energetic  components 
include:  (1)  the  radiated  energy  in  the  Geostationary  Operational  Environmental  Satellite  1-8  A  band,  (2)  the  total 
energy  radiated  from  the  soft  X-ray  (SXR)  emitting  plasma,  (3)  the  peak  energy  in  the  SXR-emitting  plasma,  (4)  the 
bolometric  radiated  energy  over  the  full  duration  of  the  event,  (5)  the  energy  in  flare-accelerated  electrons  above 
20  keV  and  in  flare-accelerated  ions  above  1  MeV,  (6)  the  kinetic  and  potential  energies  of  the  coronal  mass  ejection 
(CME),  (7)  the  energy  in  solar  energetic  particles  (SEPs)  observed  in  interplanetary  space,  and  (8)  the  amount  of 
free  (non-potential)  magnetic  energy  estimated  to  be  available  in  the  pertinent  active  region.  Major  conclusions 
include:  (1)  the  energy  radiated  by  the  SXR-emitting  plasma  exceeds,  by  about  half  an  order  of  magnitude,  the 
peak  energy  content  of  the  thermal  plasma  that  produces  this  radiation;  (2)  the  energy  content  in  flare-accelerated 
electrons  and  ions  is  sufficient  to  supply  the  bolometric  energy  radiated  across  all  wavelengths  throughout  the  event; 
(3)  the  energy  contents  of  flare-accelerated  electrons  and  ions  are  comparable;  (4)  the  energy  in  SEPs  is  typically 
a  few  percent  of  the  CME  kinetic  energy  (measured  in  the  rest  frame  of  the  solar  wind);  and  (5)  the  available 
magnetic  energy  is  sufficient  to  power  the  CME,  the  flare-accelerated  particles,  and  the  hot  thermal  plasma. 
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Sun:  X-rays,  gamma  rays 


1.  INTRODUCTION 

Solar  eruptive  events  (SEEs),  which  are  comprised  of  flares 
and  associated  coronal  mass  ejections  (CMEs),  are  the  most 
energetic  occurrences  in  the  solar  system.  Over  a  period  of  tens 
of  seconds  to  minutes,  they  can  convert  upward  of  1032  erg 
of  energy  carried  in  non-potential,  current-carrying  magnetic 
fields  into  accelerated  particles,  heated  plasma,  and  ejected  solar 
material. 

While  the  overall  energy  involved  in  a  large  SEE  is  not  in 
serious  doubt,  its  partition  among  its  component  parts  has  so  far 
been  estimated  only  for  a  few  events.  In  this  paper,  we  provide 
the  first  statistical  analysis  of  energy  partition  throughout  the 
various  manifestations  of  an  SEE,  for  38  large  events.  We 
provide  this  information  not  only  to  establish  “typical”  ratios 
of  the  energy  in  various  components  of  the  event,  but  also  to 
provide  some  idea  of  the  range  over  which  such  ratios  extend, 
and  we  especially  point  out  events  in  which  the  strength  of  one 
component  or  another  appears  to  lie  outside  the  norm.  We  offer 
this  analysis  with  the  goal  of  providing  useful  constraints  for 
modelers  of  the  energy  release  process(es)  involved. 

This  paper  grew  out  of  the  energetics  working  group  at  the 
meeting  on  “Solar  Activity  during  the  onset  of  Solar  Cycle 
24”  held  in  Napa,  CA,  from  2008  December  8-12.  It  is  a 
continuation  of  the  work  begun  at  the  Taos  ACE /RHESSI/  WIND 
joint  workshop  in  2003  that  led  to  the  works  of  Emslie  et  al. 
(2004,  2005).  These  papers  provided  the  first  detailed  analysis 
of  most  of  the  components  of  two  well-observed  SEEs  (the 
GOES  X1.5  event  on  2002  April  21  and  the  X4.8  flare  event 


of  2002  July  23),  including  the  energies  in  thermal  plasma, 
flare-accelerated  electrons  and  ions,  associated  CME,  and  solar 
energetic  particles  (SEPs).  Emslie  et  al.  (2004)  showed  that, 
for  the  two  events  in  question,  the  energy  in  the  magnetic  field 
was  sufficient  to  power  the  thermal  soft  X-ray  (SXR)  emitting 
plasma,  the  flare-accelerated  electrons  and  protons,  and  the 
kinetic  energy  in  the  CME,  and  they  also  provided  order-of- 
magnitude  estimates  of  the  partition  of  the  energy  among  these 
components.  Subsequently,  Emslie  et  al.  (2005)  also  considered 
the  energy  in  the  optical  and  EUV  continua,  and  they  cautioned 
that,  due  to  the  transfer  of  one  energy  component  to  another 
(e.g.,  flare-accelerated  electrons  thermal  plasma  -a  SXR 
emission),  care  must  be  taken  in  summing  energetic  components 
to  arrive  at  the  total  energy  released  in  an  SEE.  The  present  paper 
is  also  motivated  by  the  work  of  Mewaldt  et  al.  (2008b),  which 
was  the  first  to  address  the  ratio  of  two  energetic  components 
(the  CME  energy  in  the  rest  frame  of  the  solar  wind  and  the 
energy  in  SEPs)  for  a  statistically  significant  number  of  well- 
observed  events. 

The  basic  objective  of  the  paper  is  to  conduct  a  statistical  study 
of  the  energy  partition  into  different  components  for  many  of 
the  larger  SEEs  observed  during  the  previous  maximum  of  solar 
activity,  particularly  during  the  period  2002  February-2006 
December,  the  first  five  years  of  observations  by  the  Ramaty 
High  Energy  Solar  Spectroscopic  Imager  ( RHESSI ;  Lin  et  al. 
2002).  The  intent  is  to  apply  previously  proven  techniques  to 
determining  the  global  energetics  of  many  more  events  than  the 
two  studied  by  Emslie  et  al.  (2004,  2005),  and,  where  possible, 
to  apply  new  techniques  to  improve  the  energy  estimates. 
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Our  energy  estimates  come  from  a  wide  variety  of  observa¬ 
tions:  CME  kinetic  and  potential  energies  from  the  Large  An¬ 
gle  and  Spectrometric  COronagraph  (LASCO;  Brueckner  et  al. 
1995)  instrument  on  the  Solar  and  Heliospheric  Observatory 
(SOHO);  energy  in  flare-accelerated  charged  particles  inferred 
from  the  hard  X-rays  and  gamma  rays  observed  by  RHESSI; 
energy  contained  in  the  SXR-emitting  hot  plasma  from  the 
Geostationary  Operational  Environmental  Satellites  (GOESs) 
and  RHESSI;  energy  in  SEPs  from  the  suite  of  instruments  on 
the  Advanced  Composition  Explorer  (ACE)  and  from  GOES, 
SOHO,  the  Solar  Anomalous  and  Magnetospheric  Particle 
EXplorer  (SAMPEX),  and  the  Solar  TErrestrial  RElations 
Observatory  (STEREO);  and  total  radiated  energy  from  the 
Total  Irradiance  Monitor  (TIM;  Kopp  &  Lawrence  2005)  on  the 
SOlar  Radiation  and  Climate  Experiment  (SORCE).  For  weaker 
events,  or  where  total  irradiance  measurements  are  not  available, 
the  Flare  Irradiance  Spectral  Model  (FISM;  Chamberlin  et  al. 
2007,  2008)  was  used  to  provide  estimates  of  the  bolometric 
output  of  a  flare  based  on  other  measurements. 

In  Section  2,  we  present  the  events  studied  and  review  the 
techniques  used  to  estimate  the  different  component  energies  of 
each  event.  In  Section  3,  we  present  a  series  of  scatter  plots  of 
one  energy  component  against  another.  While  the  uncertainties 
on  the  individual  energy  estimates  are  typically  large  (often  an 
order  of  magnitude  or  greater),  these  scatter  plots,  because  of 
the  relatively  large  number  of  events  they  contain,  nevertheless 
allow  some  general  conclusions  to  be  reached  (Section  4)  about 
how  the  energy  is  partitioned.  The  spread  in  the  values  for 
the  different  energy  components  also  gives  an  idea  of  the 
uncertainties  in  the  measured  parameters  and  the  range  of  flare 
intensities  of  the  selected  events.  These  plots  also  allow  for 
the  identification  of  a  few  “outlier”  events  (Section  5)  that 
indicate  either  larger  measurement  uncertainties  or  distinctly 
different  energy  partitioning  for  those  events.  We  summarize 
the  results  in  Section  6,  which  also  provides  suggestions  for 
future  work. 

2.  COMPONENT  ENERGIES  OF  THE  SOLAR 
ERUPTIVE  EVENTS 

The  events  studied  are  listed  in  Table  1.  They  include  the 
largest  SEP  events  observed  after  2002  February,  when  RHESSI 
was  launched,  excluding  those  events  beyond  the  west  limb  (for 
which  no  reliable  active  region  identification  can  be  made)  and 
those  events  located  from  E60°  to  E90°  (for  which  the  evaluation 
of  the  SEP  energy  is  highly  uncertain — see  Section  2.7).  They 
also  include  the  two  events  studied  by  Emslie  et  al.  (2004), 
which  appear  as  Events  2  and  6  in  Table  1.  Additional  events 
include  all  flares  for  which  RHESSI  detected  significant  (>4er) 
emission  in  the  2.223  MeV  neutron-capture  gamma-ray  line 
(Shih  2009;  Shih  et  al.  2009).  This,  plus  the  inclusion  of  an 
intriguing  behind-the-limb  event  with  a  strong  CME  on  2002 
July  20  (Event  5  in  Table  1),  resulted  in  a  total  of  38  events  for 
study.  As  permitted  by  the  available  data,  estimates  were  made 
of  the  following  energies  for  each  of  the  38  events. 

1.  Radiated  energy  in  the  GOES  1-8  A  band. 

2.  Total  radiated  energy  from  the  SXR-emitting  plasma. 

3.  Total  (bolometric)  radiated  output. 

4.  Peak  thermal  energy  of  the  SXR-emitting  plasma. 

5.  Energy  in  flare-accelerated  electrons. 

6.  Energy  in  flare-accelerated  ions. 

7.  CME  kinetic  energy  in  the  rest  frame  of  the  Sun. 

8.  CME  kinetic  energy  in  the  solar  wind  rest  frame. 


9.  CME  gravitational  potential  energy. 

10.  Energy  in  SEPs. 

11.  Free  (non-potential)  magnetic  energy  in  the  active  region. 

It  is  important  to  keep  in  mind  the  differences  among  the 
first  four  items  on  this  list.  They  are  all  related,  but  are 
included  separately  since  they  can  each  be  estimated  indepen¬ 
dently,  and  indeed  relatively  straightforwardly,  from  the  avail¬ 
able  measurements,  and  since  collectively  they  provide  signif¬ 
icant  information  on  the  thermal  energy  of  each  flare,  how  it 
is  distributed  in  temperature,  and  when  it  is  generated  and 
released.  Further  details  on  these  four  items  are  provided  in 
Sections  2. 1-2.3.  Broadly  speaking,  the  first  item  is  the  energy 
radiated  in  the  narrow  GOES  band  from  1  to  8  A,  obtained 
directly  from  background-subtracted  data  (Section  2.1).  The 
second  is  the  energy  radiated  over  all  wavelengths  (including 
the  1-8  A  band)  from  the  hot  SXR-emitting  plasma  and  is  a 
quantity  inferred  from  the  plasma  parameters  (emission  measure 
and  temperature)  revealed  by  the  GOES  1-8  A  measurements. 
The  third  item  is  the  total  energy  radiated  over  all  wavelengths 
from  all  components  of  the  flare  at  all  temperatures  (including 
that  from  the  SXR-emitting  plasma);  in  some  cases  this  is  di¬ 
rectly  observed  and  in  some  cases  inferred  from  modeling  of 
the  emission  in  select  wavelength  ranges — see  Section  2.2. 
The  fourth  item  does  not  specify  a  radiated  energy  at 
all,  but  rather  the  peak  thermal  energy  content  of  the  hot 
SXR-emitting  plasma;  this  quantity  is  inferred  (Section  2.3) 
from  the  parameters  of  spectroscopic  fits  to  RHESSI  data.  It  is 
important  to  realize  (Emslie  et  al.  2005)  that  these  four  com¬ 
ponents  are  not  separate  flare  energy  components  and  therefore 
cannot  be  summed  together  to  obtain  a  total  flare  thermal  energy. 

We  have  not  evaluated  energy  losses  from  the  SXR-emitting 
plasma  by  thermal  conduction.  However,  we  would  note  that 
conductive  transfer  of  thermal  energy  from  hot  SXR-emitting 
plasma  into  the  relatively  cool  chromospheric  plasma  will 
generally  result  in  the  thermal  energy  content  of  the  hot 
SXR-emitting  plasma  being  quite  effectively  transported  to,  and 
ultimately  radiated  away  by,  such  relatively  cool  plasma,  one 
contribution  to  the  total  (bolometric)  radiated  output — the  third 
item  on  the  list.  Further,  for  the  events  considered  here  there 
is  little  observational  data  available  on  the  energy  contained 
in  turbulence  and  directed  mass  motions  of  thermal  plasma, 
components  that  may  well  contain  energies  comparable  to  the 
thermal  energy  of  the  SXR-emitting  plasma  (see,  e.g.,  Doschek 
et  al.  1992). 

The  data  on  all  the  component  energies  are  summarized 
in  Table  1.  GOES  SXR  data  are  available  for  all  the  events. 
However,  because  of  missing  or  inadequate  data  or  limited 
instrument  sensitivities,  definitive  energy  estimates  for  all  the 
energy  components  listed  above  are  available  for  only  six  events 
(Events  13,  14,  20,  23,  25,  and  38).  As  mentioned  above,  Event 
5  was  located  behind  the  limb;  thus  only  a  “plausible”  magnetic 
energy  estimate  (not  included  in  Table  1 )  could  be  obtained  from 
observations  of  the  most  likely  responsible  active  region  once  it 
had  moved  onto  the  solar  disk,  and  the  listed  radiated  energies 
are  lower  limits. 

2.1.  Radiated  Energy  from  Hot  Plasma 

For  each  event,  we  estimated  the  time-integrated  SXR  and 
total  radiated  energies  from  the  hot  SXR-emitting  plasma — the 
columns  labeled  “SXR”  and  “T-rad,”  respectively,  in  Table  1. 
Fluxes  in  W  m 1  2 3 4 5 6 7 8  for  the  1-8  A  and  0.5-4  A  bands  are  provided 
by  one  of  NOAA  GOES  satellites  every  3  s.  The  total  emission 
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Table  1 

Event  List  with  Component  Energies  ( x  1 0 10  erg) 


No. 

Datea 

Timeb 

Class 

SXRC 

T-radd 

BoP 

Peakf 

Elec® 

Ionh 

KE* 1 

SWJ 

PEk 

SEP1 

Magm 

1 

02/02/20 

05:52 

M5.1 

0.043 

1.2 

13 

17 

5.6 

6.3 

0.13 

1200 

2 

02/04/21 

00:43 

XI. 5 

1.2 

38 

150 

13 

20 

230 

160 

5.0 

23 

660 

3 

02/05/22 

03:18 

C5.0 

0.048 

5.6 

9 

84 

45 

10 

2.7 

260 

4 

02/07/15 

19:59 

X3.0 

0.31 

6.4 

44 

>2.2 

>3.6 

160 

76 

10 

3.8 

1500 

5 

02/07/20 

21:04 

X3.3n 

>1.5 

>26 

>210 

260 

170 

6 

02/07/23 

00:18 

X4.8 

1.2 

19 

150 

2.5 

32 

39 

260 

150 

20 

<30 

2000 

7 

02/08/24 

00:49 

X3.1 

1.1 

24 

160 

5.9 

11 

210 

130 

16 

3.9 

2500 

8 

02/11/09 

13:08 

M4.6 

0.11 

5.0 

8 

1.3 

60 

180 

110 

20 

0.51 

550 

9 

03/05/27 

22:56 

X1.4 

0.16 

3.6 

16 

2.8 

7.4 

0.19 

260 

10 

03/06/17 

22:27 

M6.9 

0.21 

4.6 

17 

2.4 

4.6 

6.7 

140 

11 

03/10/26 

17:21 

XI. 2 

1.2 

31 

88 

240 

130 

32 

0.75 

1700 

12 

03/10/28 

09:51 

X17 

4.4 

68 

362° 

>19 

>56 

>190 

1200 

850 

63 

43 

2900 

13 

03/10/29 

20:37 

X10 

1.9 

31 

137° 

11 

no 

30 

340 

220 

25 

9.7 

2900 

14 

03/11/02 

17:03 

X8.3 

1.8 

24 

130 

9.3 

130 

68 

270 

200 

10 

9.3 

2800 

15 

03/11/03 

09:43 

X3.9 

1.1 

17 

97 

2.4 

120 

3.1 

780 

16 

03/11/04 

19:29 

X28 

4.8 

72 

426° 

>3.1 

>21 

610 

410 

25 

5.3 

2800 

17 

04/07/15 

18:15 

XI. 6 

0.16 

4.1 

8 

0.93 

42 

<0.1 

820 

18 

04/07/25 

05:39 

M7.1 

0.069 

1.3 

10 

2.9 

2300 

19 

04/11/07 

15:42 

X2.0 

0.32 

5.0 

56 

3.0 

43 

220 

130 

25 

4.2 

610 

20 

04/11/10 

01:59 

X2.5 

0.32 

7.7 

15 

2.0 

20 

3.4 

230 

180 

16 

2.4 

610 

21 

05/01/15 

00:22 

XI. 2 

0.23 

4.7 

23 

5.0 

32 

1500 

22 

05/01/15 

22:25 

X2.6 

1.3 

22 

78 

7.1 

63 

15 

730 

540 

1600 

23 

05/01/17 

06:59 

X3.8 

1.8 

34 

150 

17 

48 

13 

1000 

730 

50 

11 

1600 

24 

05/01/19 

08:03 

XI. 3 

0.43 

7.0 

54 

5.9 

82 

29 

1600 

25 

05/01/20 

06:36 

X7.1 

2.9 

43 

150 

10 

25 

120 

15-79 

7.8-61 

2.0 

7.8 

1600 

26 

05/05/13 

16:13 

M8.0 

0.44 

14 

49 

3.1 

13 

39 

22 

4.0 

7.3 

400 

27 

05/07/14 

10:16 

XI. 2 

0.64 

12 

87 

4.3 

24 

100 

66 

6.3 

2.9 

310 

28 

05/07/27 

04:33 

M3. 7 

0.16 

4.5 

30 

1.3 

12 

100 

62 

10 

310 

29 

05/08/22 

16:46 

M5.6 

0.34 

9.8 

35 

3.2 

6.3 

110 

76 

10 

6.4 

390 

30 

05/08/25 

04:31 

M6.4 

0.050 

1.2 

11 

1.1 

16 

<1.9 

110 

31 

05/09/07 

17:17 

X17 

4.9 

68 

322° 

>5.6 

>10 

>0.7 

1400 

32 

05/09/09 

19:13 

X6.2 

3.1 

44 

250 

>7.9 

>120 

>1.7 

1300 

33 

05/09/10 

21:30 

X2.1 

0.99 

17 

82 

6.0 

13 

1.0 

1300 

34 

05/09/13 

19:19 

X1.5 

1.1 

25 

85 

330 

200 

32 

3.0 

1400 

35 

05/09/13 

23:15 

XI. 7 

0.23 

4.7 

21 

2.3 

32 

<0.1 

1400 

36 

06/12/05 

10:18 

X9.0 

1.4 

19 

92 

>5.1 

>360 

>4.5 

400 

37 

06/12/06 

18:29 

X6.5 

1.1 

18 

59° 

6.8 

40 

36 

410 

38 

06/12/13 

02:14 

X3.0 

1.1 

17 

75 

4.8 

13 

14 

74 

44 

6.3 

3.2 

570 

Notes. 

a  In  yy/mm/dd  format. 

b  GOES  start  time  (UT). 

c  Radiated  energy  in  the  GOES  1-8  A  band. 

d  Total  radiated  energy  from  the  SXR-emitting  plasma. 

e  Bolometric  radiated  energy. 

f  Peak  thermal  energy  of  the  SXR-emitting  plasma. 

8  Energy  in  flare-accelerated  electrons. 
h  Energy  in  flare-accelerated  ions. 

1  CME  kinetic  energy  in  the  rest  frame  of  the  Sun. 

1  CME  kinetic  energy  in  the  solar  wind  rest  frame. 
k  CME  gravitational  potential  energy. 

1  Energy  in  SEPs. 

m  Non-potential  magnetic  energy  in  the  active  region. 
n  Behind-the-limb  event. 

0  Bolometric  irradiance  directly  measured  with  TIM — see  Table  2. 


in  the  1-8  A  band  (“SXR”)  is  obtained  simply  by  summing 
the  background-subtracted  fluxes  over  the  duration  of  the  flare, 
from  the  GOES  start  time  (given  by  NOAA  and  listed  in  Table  1) 
to  the  time  when  the  flux  had  decreased  to  10%  of  the  peak  value. 
The  background  that  was  subtracted  was  taken  as  the  lowest  flux 
in  the  hour  or  so  before  and/or  after  the  flare. 

To  calculate  the  radiated  energy  from  the  hot  plasma,  we 
used  the  measured  GOES  SXR  fluxes  in  a  manner  similar  to 


that  described  in  Emslie  et  al.  (2004),  specifically  using  the 
1DL  GOES  Workbench  available  in  SolarSoftware  (SSW).  This 
allows  us  to  obtain  a  consistent  set  of  values  for  all  events 
since  GOES,  unlike  RIIESSI,  has  full  coverage  for  all  events. 
This  calculation  assumes  that  the  hot  plasma  at  any  given 
time  is  isothermal;  the  temperature  and  emission  measure  are 
calculated  from  the  two-channel  GOES  data  using  the  relations 
given  by  White  et  al.  (2005).  Using  the  emission  measure 


3 


The  Astrophysical  Journal,  759:71  (18pp),  2012  November  1 


Emslie  et  al. 


Table  2 

TIM  and  FISM  Bolometric  Energies  (1030  erg) 


Event 

No. 

Date 

Total3 

Uncertainty3 

TIM 

Revised  Estimate 

Corrected15 

FISM 

Corrected15 

Difference 

(TIM  -  FISM) /TIM 

12 

2003  Oct  28 

600 

39% 

362 

362 

310 

14% 

13 

2003  Oct  29 

240 

86% 

137 

137 

128 

7% 

16 

2003  Nov  4 

260 

65% 

142 

426 

447 

-5% 

31 

2005  Sep  7 

300 

71% 

150 

322 

266 

17% 

37 

2006  Dec  6 

46c 

59 

82 

-39% 

Notes. 

a  Woods  et  al.  (2006). 

b  Corrected  for  limb  darkening. 

c  C.  S.  Moore  et  al.  (2012,  in  preparation). 

and  temperature  thus  inferred,  and  the  optically  thin  radiation 
loss  rate  versus  temperature  function  (for  coronal  abundances 
and  Mazzotta  et  al.  1998  ionization  equilibria)  taken  from  the 
CHIANTI  database  (Dere  et  al.  1997,  2009),  we  used  the  IDL 
procedure  rad_loss,  available  in  SSW,  to  calculate,  for  each 
3  s  time  interval,  the  energy  radiated  from  the  SXR-emitting 
plasma  over  all  wavelengths.  Finally,  we  summed  the  radiated 
energies  over  the  duration  of  the  flare  (from  the  GOES  start 
time  until  the  1-8  A  flux  decreased  to  10%  of  its  peak  value)  to 
produce  the  total  radiated  energy  given  in  the  column  labeled 
“T-rad”  in  Table  1 .  Significant  energy  could  be  radiated  after  this 
nominal  end  of  the  flare,  particularly  if  there  is  a  “second  phase” 
that,  according  to  Woods  et  al.  (2011)  and  Su  et  al.  (2012),  can 
release  an  amount  of  energy  that  is  similar  to  that  released  in  the 
initial  phase.  Generally,  however,  the  values  quoted  in  Table  1 
should  include  more  than  50%  of  the  energy  radiated  by  the 
SXR-emitting  plasma. 

2.2.  Bolometric  Irradiance 

Estimates  of  the  bolometric  irradiance,  the  total  energy  radi¬ 
ated  from  the  flare  integrated  across  the  entire  solar  spectrum, 
for  each  of  the  events  are  provided  in  the  column  labeled  “Bol” 
in  Table  1 .  For  five  of  the  events  listed  in  Table  1 ,  the  bolometric 
irradiance  was  measured  directly  by  the  TIM  (Kopp  &  Lawrence 
2005)  on  board  SORCE  as  an  increase  in  the  total  solar  irradi¬ 
ance  above  the  (highly  variable)  pre-  and  post-flare  background 
levels.  Total  flare  irradiance  values  were  reported  in  this  manner 
for  Events  12,  13,  16,  and  31  by  Woods  et  al.  (2006),  and  the 
bolometric  irradiance  for  Event  37  will  be  reported  by  C.  S. 
Moore  et  al.  (2012,  in  preparation).  Both  previously  published 
values  and  revised  estimates  made  for  this  paper  are  listed  in 
Table  2.  A  final  correction  factor  (see  Table  2)  was  then  applied 
to  allow  for  limb-darkening  absorption  when  the  path  to  the  ob¬ 
server  becomes  optically  thick  at  some  wavelengths;  the  value 
of  this  factor  can  be  up  to  ~3.0 — see  Equation  (2)  in  Woods 
et  al.  (2006). 

To  complement  these  direct  measurements  and  so  provide 
a  consistent  set  of  bolometric  values  for  all  of  the  events  in 
Table  1,  estimates  from  the  FISM  (Chamberlin  et  al.  2007, 
2008)  were  used,  with  various  assumptions  and  corrections  as 
described  below.  FISM  is  an  empirical  model  that  provides 
estimates  of  the  total  amount  of  solar  radiated  energy  over  a 
broad  wavelength  range  from  1  to  1900  A  and  over  a  wide 
range  of  timescales  from  seconds  to  years.  It  uses  measure¬ 
ments  in  this  wavelength  range  from  the  Solar  EUV  Experiment 
(SEE;  Woods  et  al.  2005)  on  the  Thermosphere  Ionosphere 
Mesosphere  Energetics  and  Dynamics  {TIMED)  satellite 
and  the  SOLar-STellar  Irradiance  Comparison  Experiment 


(SOLSTICE;  Rottman  et  al.  1993)  on  the  Upper  Atmosphere 
Research  Satellite. 

For  the  relatively  rapid  GOES  1-8  A  SXR  flux  variations  that 
occur  during  a  solar  flare,  different  empirical  factors  appropriate 
to  the  rise  and  decay  phases  of  the  flare,  respectively,  are  used 
to  relate  the  SXR  flux  to  the  total  radiated  energy  during  those 
phases.  Various  daily  proxies  are  also  used  to  represent  the 
more  gradual  variations  in  solar  irradiance  due  to  active  region 
evolution,  solar  rotation,  and  the  solar  cycle.  The  daily  pre-flare 
irradiance  spectrum  is  subtracted  from  each  value  to  get  the 
radiated  energy  from  the  flare  alone,  and  this  is  then  integrated 
over  the  duration  of  the  GOES  flare  to  get  the  total  radiated 
energy  in  erg  cm-2  at  the  detector.  Then,  assuming  uniform 
radiation  over  2n  steradians,  the  total  radiated  energy  from  the 
flare  in  the  1-1900  A  wavelength  range  can  be  calculated,  with 
1  minute  cadence. 

The  1-1900  A  solar  irradiance  is  converted  to  total  radiated 
energy  over  all  wavelengths  (the  bolometric  irradiance)  by 
multiplying  by  an  empirical  conversion  factor  of  2.42  ±  0.31, 
determined  by  comparing  the  1-1900  A  solar  irradiance  with 
the  absolute  bolometric  intensity  for  the  five  flares  for  which 
the  latter  could  be  measured  directly  (see  Table  2). 

The  uncertainties  on  the  calculated  values  of  the  bolometric 
irradiance  listed  in  Tables  1  and  2  are  made  up  of  several  parts. 
The  dominant  uncertainty  comes  from  the  TIM  measurements 
themselves  and  is  due  to  the  variations  in  the  total  solar 
irradiance  of  the  non-flaring  Sun.  Other  contributions  to  the 
overall  uncertainty  are  the  errors  on  the  FISM  estimates  of 
the  1-1900  A  flux,  the  conversion  from  UV  irradiance  to  total 
solar  irradiance,  and  the  limb-darkening  correction.  The  overall 
uncertainty  on  the  calculated  values  is  ±  ~70%  for  those  events 
that  are  near  disk  center  and  ±  ~90%  for  the  near-limb  events.8 

2.3.  Peak  Thermal  Energy  Content  of  the  Hot  Plasma 

The  peak  thermal  energy  content  of  the  SXR-emitting  plasma 
(the  column  labeled  “Peak”  in  Table  1)  was  determined  from 
RHESSI  imaging  spectroscopy  data.  We  first  fit  the  observed 
RHESSI  hard  X-ray  spectra  with  the  sum  of  a  single-temperature 
Maxwellian  plus  the  form  expected  from  a  double-power-law 
electron  spectrum  (Equation  (2)  in  Section  2.4).  The  fit  param¬ 
eters  appropriate  to  both  thermal  and  nonthermal  components 
were  determined  for  each  time  interval  using  the  forward-fitting 
method  implemented  in  the  OSPEX  software  package  available 
in  SSW.  The  temperatures  and  emission  measures  obtained  from 


8  Because  the  conversion  of  the  FISM  radiated  energy  to  bolometric  energy  is 
based  on  the  five  events  measured  directly  with  TIM,  the  bolometric  energies 
for  these  five  events  derived  from  the  FISM  estimate  differ  (after  correcting  for 
limb  darkening)  by  less  than  40%  from  the  directly  measured  values. 
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RHESSI  in  this  way  tend  to  agree  closely  with  the  corresponding 
values  obtained  from  the  standard  GOES  data  analysis  discussed 
in  Section  2.1.  However,  somewhat  higher  temperatures  can  be 
obtained  because  of  the  RHESSI  coverage  to  higher  energies;  in¬ 
deed,  superhot  components  with  reported  temperatures  as  high 
as  ~50  MK  may  exist  in  some  flares  and  are  not  accurately 
reflected  in  the  GOES  thermal  analysis  (Lin  et  al.  1981;  Caspi 
2010;  Caspi  &  Lin  2010).  The  thermal  function  included  the 
line-plus-continuum  components  determined  using  CHIANTI, 
again  with  coronal  abundances  and  Mazzotta  et  al.  (1998)  ion¬ 
ization  equilibria.  From  these  fits,  the  average  temperature  Tq 
(K)  and  emission  measure  EM  =  /  n2  dV  (cm-3)  of  the  thermal 
plasma  were  determined  every  20  s  throughout  the  flare.  (Here 
ne  is  the  electron  density  (cm-3)  and  L  is  the  emitting  volume 
(cm3).) 

The  thermal  energy  content  f/th  of  the  plasma  can  be  calcu¬ 
lated  from  the  expression  (e.g.,  de  Jager  et  al.  1986) 

£4  =  3  ne  kT0  f  Lap  ~  4.14  x  10-16  T0  ^EM  x  /  Lap  erg, 

(1) 


where  k  is  Boltzmann’s  constant,  /  is  the  volumetric  filling  fac¬ 
tor,  and  14p  is  the  apparent  volume  of  the  SXR  source.  Starting 
in  2003,  SXR  images  are  also  available  from  the  GOES  Soft 
X-ray  Imagers  (SXIs;  Lemen  et  al.  2004;  Pizzo  et  al.  2005); 
such  images  could  be  used  to  provide  estimates  of  Lap,  as  de¬ 
scribed  in  Emslie  et  al.  (2004).  However,  both  for  consistency 
with  the  earlier  analysis  of  Emslie  et  al.  (2004,  which  analyzed 
events  that  occurred  in  2002,  prior  to  the  SXI  deployment), 
and  since  the  parameters  EM  and  Tq  in  Equation  (1)  are  de¬ 
duced  from  RHESSI  data,  we  have  chosen  to  use  estimates  of 
Lap  that  are  deduced  from  RHESSI  images,  made  using  the  3a 
clean  method  of  Dennis  &  Pernak  (2009).  We  further  take /to 
be  unity,  consistent  with  Emslie  et  al.  (2004)  and  further  jus¬ 
tified  by  the  recent  work  of  Guo  et  al.  (2012),  who  used  hard 
X-ray  imaging  spectroscopy  data  of  22  extended-loop  events  to 
derive  a  (logarithmic)  mean  filling  factor  /  =  0.20  x/-t-  3.9  (la 
standard  error). 

The  peak  energy  values  listed  in  Table  1  are  the  highest  values 
of  Uth  obtained  from  this  analysis,  usually  at  or  near  the  time  of 
the  peak  GOES  flux. 

2.4.  Flare-accelerated  Electrons 


The  energies  in  flare-accelerated  electrons  are  listed  in  the 
column  labeled  “Elec”  in  Table  1.  They  were  determined  by 
using  the  OSPEX  algorithm  to  fit  a  combined  isothermal-plus- 
nonthermal  function  to  the  measured  RHESSI  spatially  inte¬ 
grated  X-ray  spectra.  The  nonthermal  component  was  assumed 
to  be  bremsstrahlung  from  energetic  electrons  with  an  injected 
spectrum  Fq(Eq)  (electrons  cm-2  s-1  keV-1)  in  the  form  of  a 
broken  power  law: 


F0(E0)  =  A 


0, 

(e0/ep)-s\ 

(E0/EprsHEb/Ep)s^, 

0, 


Eq  <  Emin 
Emin  ^  Eo  <  Eb 
l-h  L  Eo  <  Emax 
Emax  ^  E 0- 

(2) 


The  seven  parameters  of  this  model  spectrum  are  the  normal¬ 
ization  parameter  A,  the  low-  and  high-energy  cutoffs  Emin  and 
£max,  the  break  energy  Eb,  and  the  power-law  indices  ilq  and 
82  below  and  above  the  break  energy,  respectively.  The  (arbi¬ 
trary)  value  of  the  pivot  energy  Ep  was  fixed  at  50  keV.  Also, 
the  high-energy  cutoff  Emax  was  fixed  at  30  MeV,  an  energy  so 


high  above  the  energy  range  of  interest  (<500  keV)  that  it  has 
a  negligible  effect  on  the  calculated  X-ray  spectrum  and  so  is 
equivalent  to  having  no  high-energy  cutoff  at  all. 

The  OSPEX  analysis  uses  a  forward-fitting  procedure  that 
starts  by  dividing  the  flare  into  multiple  time  intervals — here 
we  used  20  s  intervals.  For  each  interval,  the  function  thick2  in 
SSW  is  used  with  a  set  of  starting  parameters  for  the  electron 
spectrum  (2)  used  to  calculate  the  X-ray  photon  spectrum, 
assuming  electron-ion  bremsstrahlung  in  a  thick  target  that  is 
“cold”  in  the  sense  that  the  ambient  electrons  have  a  mean 
energy  kT  that  is  significantly  lower  than  the  lowest  energy 
of  the  accelerated  electrons.  In  general,  consideration  must 
also  be  given  to  the  ionization  state  of  the  target,  since  the 
bremsstrahlung  efficiency  is  a  factor  of  ~3  times  higher  for  a 
fully  ionized  plasma  than  for  an  un-ionized  gas  (Brown  1973; 
Kontar  et  al.  2003).  However,  since  most  of  the  beam  energy 
is  in  the  lower  energy  electrons  that  stop  higher  in  the  corona, 
we  used  parameters  appropriate  for  a  fully  ionized  plasma  to 
calculate  the  total  nonthermal  energy.  A  more  refined  calculation 
is  possible  using  the  procedure  outlined  by  Kontar  et  al.  (2002) 
and  Su  et  al.  (201 1),  but  no  significant  difference  is  expected  in 
the  resulting  total  energy  in  electrons  above  Emin- 

The  resulting  photon  spectrum  is  then  folded  through  the 
detector  response  matrix  to  generate  a  count-rate  spectrum, 
which  is  added  to  the  count-rate  spectrum  calculated  for  the 
thermal  spectrum  discussed  in  Section  2.1.  Then,  through  an 
iterative  procedure,  we  find  best-fit  values  of  the  parameters 
describing  the  electron  spectrum  (2)  by  minimizing  the  y2 
statistic  between  the  calculated  and  the  measured  background- 
subtracted  count  spectra.  The  total  energy  Ue  in  electrons  for  a 
given  event  is  then  computed  by  integrating  the  best-fit  electron 
energy  spectrum  above  £min  for  each  time  interval  and  summing 
the  results  over  the  duration  of  the  flare,  resulting  in  the  values 
listed  in  the  column  labeled  “Elec”  in  Table  1. 

In  order  to  obtain  the  most  reliable  spectral  fits  to  the  RHESSI 
data  and  thus  better  evaluate  the  uncertainties  in  the  calculated 
values  of  Ue ,  we  chose  to  use  data  from  just  one  of  RHESSI’ s  nine 
detectors — detector  4.  This  particular  detector  has  good  energy 
resolution  and  sensitivity,  which  allowed  us  to  apply  the  most 
up-to-date  corrections  for  energy  resolution  and  calibration, 
photospheric  albedo,  pulse  pile-up,  and  background  subtraction 
that  are  available  with  the  current  analysis  software.  For  the 
large  events  studied,  the  count  rates  were  sufficiently  high  that 
selecting  just  a  single  detector  did  not  seriously  degrade  the 
spectroscopic  capability  up  to  the  photon  energies  required  to 
determine  the  parameters  of  interest.  Milligan  &  Dennis  (2009) 
have  shown  that  similar  best-fit  parameter  values  are  determined 
using  different  individual  detectors  (detectors  1,  3,  4,  5,  6,  and 
9  in  their  case),  which  leads  to  an  estimate  of  the  systematic 
uncertainties  in  the  calculated  total  energy  in  electrons  of  ~20%. 
This  is  negligible  compared  to  the  uncertainty  arising  from  the 
difficulty  in  establishing  the  value  of  the  low-energy  cutoff 
energy  Enlln,  an  uncertainty  that  arises  because  the  thermal 
emission  generally  dominates  the  low-energy  part  of  the  X-ray 
photon  spectrum  up  to  energies  where  the  effects  of  a  cutoff  in 
the  electron  spectrum  might  be  detectable.  We  used  the  largest 
value  of  £min  that  still  gave  an  acceptable  fit  (reduced  y1  ~  1). 
As  a  result,  the  values  of  Ue  listed  in  Table  1  are  lower  limits 
to  the  energy  in  the  nonthermal  electrons.  Furthermore,  because 
of  the  steep  form  of  the  electron  spectra  (<5j  >  4),  these  values 
are  particularly  sensitive  to  £mjn,  so  that  the  energies  in  flare- 
accelerated  electrons  could  be  up  to  an  order  of  magnitude 
higher  than  those  reported  in  Table  1 . 
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2.5.  Flare-accelerated  Ions 

The  energies  in  flare-accelerated  ions  with  energies  above 
1  MeV  are  listed  in  the  column  labeled  “Ion”  in  Table  1. 
In  order  to  provide  a  consistent  set  of  values  for  as  many 
events  as  possible,  the  energies  were  estimated  solely  from 
RHESSI  measurements  of  the  fluence  (time  integral  of  the  flux, 
photons  cm-2)  in  the  2.223  MeV  neutron-capture  gamma-ray 
line.  Our  sample  of  events  is  primarily  based  on  the  studies  of 
Shih  (2009)  and  Shih  et  al.  (2009),  who  analyzed  RHESSI  flares 
from  2002  to  2005  that  had  either  2.223  MeV  line  emission 
and/or  >0.3  MeV  electron  bremsstrahlung  continuum  emis¬ 
sion.  Of  those  flares,  energies  are  included  only  for  those  that 
have  >2(7  detections  of  that  line,  and  further  only  as  4 a  upper 
limits  if  below  a  4 a  detection.  We  also  include  three  additional 
flares  that  occurred  in  2006  (Events  36-38).  We  chose  a  lower 
energy  threshold  of  1  MeV  because  the  production  of  detected 
nuclear  gamma-ray  lines  from  elements  such  as  20Ne  begins  at 
energies  as  low  as  ~3  MeV,  and  it  is  therefore  evident  that  the 
ion  spectrum  extends  down  to  ~1  MeV,  at  least  in  a  few  large 
events  (Ramaty  et  al.  1995;  Ramaty  &  Mandzhavidze  2000). 
The  spectral  shape  is  essentially  unknown  below  1  MeV. 

In  order  to  estimate  the  energy  in  ions  from  the  2.223  MeV 
line  fluences,  the  following  steps  were  taken.  The  measured 
fluences  of  the  line  were  first  corrected  for  attenuation  in  the 
solar  atmosphere  assuming  a  given  depth  of  production  of  the 
photons  (Hua  &  Lingenfelter  1987)  and  allowing  for  the  flare 
position  on  the  solar  disk.  The  corrected  fluence  values  were  then 
converted  to  the  proton  energy  above  30  MeV  using  conversion 
factors  given  by  Murphy  et  al.  (2007)  and  Shih  (2009).  The 
30  MeV  threshold  was  used  at  this  stage  in  the  analysis  because 
the  2.223  MeV  line  is  produced  by  ions  with  energies  >20  MeV 
nucleon-1,  so  that  the  conversion  factors  are  less  dependent  on 
the  assumed  power-law  index  of  the  proton  spectrum. 

In  order  to  estimate  the  energy  in  protons  above  1  MeV,  an 
extrapolation  is  required  over  one-and-a-half  orders  of  magni¬ 
tude  in  proton  energy,  so  that  the  inferred  energy  above  1  MeV 
depends  critically  on  the  spectral  index  used  in  this  extrapo¬ 
lation.  For  the  largest  RHESSI  flares,  where  multiple  types  of 
ion-associated  gamma-ray  emission  can  be  detected  and  fitted 
simultaneously,  the  ion  power-law  spectral  indices  are  found 
to  be  typically  in  the  range  3-5,  a  range  of  indices  consistent 
with  that  found  in  a  study  of  flares  observed  by  the  Gamma-Ray 
Spectrometer  on  the  Solar  Maximum  Mission  (Ramaty  et  al. 
1996).  Consequently,  for  the  purposes  of  estimating  the  total 
energy  in  protons,  we  have  assumed  a  power-law  proton  spec¬ 
trum  with  a  single  spectral  index  of  4  that  extends  down  to  a 
lower  cutoff  energy  of  1  MeV.  Because  of  the  long  “lever  arm” 
associated  with  this  extrapolation,  an  uncertainty  in  the  spectral 
index  of  ±1  corresponds  to  an  uncertainty  in  the  total  energy 
content  above  1  MeV  of  about  ±1.5  orders  of  magnitude. 

Even  under  the  assumption  that  the  spectra  for  the  various 
types  of  ions  have  the  same  spectral  index  and  low-energy  cutoff, 
the  energy  content  will  also  depend  on  the  accelerated  particle 
composition.  The  ratio  of  the  energy  content  in  all  ions  (protons 
plus  a-particles  and  heavier  nuclei)  to  the  energy  in  protons 
can  vary  between  ~2  and  ~6;  here  we  assume  that  the  energy 
in  flare-accelerated  ions  is  three  times  the  energy  content  in 
flare-accelerated  protons. 

For  a  number  of  the  events,  the  total  energy  content  of  ions 
listed  in  Table  1  is  a  lower  limit  because  RHESSI  did  not  see  the 
complete  time  history  as  the  result  of  spacecraft  night  or  passage 
through  the  South  Atlantic  Anomaly.  In  addition,  there  can  be 


other  complications  that  affect  the  observation  or  interpretation 
of  the  neutron-capture  line.  The  affected  events  are  as  follows. 

1.  Event  12.  RHESSI  missed  a  significant  fraction  of  the 
neutron-capture  line  emission,  including  the  peak,  as  shown 
by  observations  of  this  flare  by  International  Gamma-Ray 
Astrophysics  Laboratory  (Kiener  et  al.  2006). 

2.  Event  31:  RHESSI  observed  only  ~2  minutes  of  a  signif¬ 
icantly  longer  impulsive  phase.  Furthermore,  the  level  of 
atmospheric  attenuation  is  very  uncertain  due  to  this  flare’s 
large  heliocentric  angle,  so  we  use  a  conservative  angle  of 
80°  to  determine  the  correction  factor. 

3.  Events  32  and  37.  RHESSI  missed  the  peak  of  the  impulsive 
emission,  and  thus  possibly  a  significant  fraction  of  the  total 
emission. 

4.  Event  36.  RHESSI  missed  some  fraction  of  the  2.223  MeV 
emission  as  it  was  just  coming  out  of  Earth  shadow.  RHESSI 
observations  started  at  10:31  UT  and  the  GOES  X-ray  flare 
started  at  ~  10: 19  UT. 

5.  Events  14,  15,  22,  33,  and  38.  RHESSI  likely  missed  a 
small  fraction  of  the  neutron-capture  line  emission  late  in 
the  flare.  However,  this  missing  energy  is  smaller  than  the 
other  uncertainties  in  the  energy  estimates  discussed  above. 

6.  Events  36-38.  By  2006  December,  RHESSF s  detectors  had 
reduced  gamma-ray  sensitivity  resulting  from  accumulated 
radiation  damage,  a  reduction  that  is  difficult  to  estimate. 

2.6.  Coronal  Mass  Ejection 

The  CME  kinetic  energies,  both  in  the  rest  frame  of  the  Sun 
and  in  that  of  the  solar  wind  (for  comparison  with  the  SEP 
energies — see  Section  3.6),  are  listed  in  the  columns  labeled 
“KE”  and  “SW,”  respectively,  of  Table  1.  The  gravitational 
potential  energies  of  the  CMEs  are  listed  in  the  column  labeled 
“PE.”  These  CME  energies  were  estimated  from  calibrated 
LASCO  images  using  the  procedure  detailed  in  Vourlidas  et  al. 
(2010,2011). 

Briefly,  this  procedure  consists  of  the  following  steps.  First, 
we  selected  two  LASCO  images,  one  containing  the  CME  and 
the  other  taken  before  the  event  as  close  in  time  as  possible 
to  the  flare  with  no  disturbances  or  ejecta  over  the  path  of  the 
subsequent  CME.  Next,  the  images  were  calibrated  (in  units  of 
mean  solar  brightness)  and  the  pre-event  image  was  subtracted 
from  the  CME  image.  The  excess  brightness  revealed  by  this 
subtracted  image  is  due  to  Thomson  scattering  of  photospheric 
radiation  from  the  excess  mass  in  the  CME.  This  excess 
brightness  can  therefore  be  converted  to  excess  mass  of  the 
CME  under  the  usual  assumptions  that  (1)  all  of  the  CME  mass 
is  concentrated  on  the  plane  of  the  sky  and  (2)  the  CME  material 
consists  of  90%  H  and  10%  He  (Poland  et  al.  1981;  Vourlidas 
et  al.  2000,  2010).  We  used  the  first  assumption  because  the  true 
three-dimensional  distribution  of  the  CME  mass  along  the  line 
of  sight  (LOS)  is  unknown.  The  second  assumption  represents 
an  “average”  coronal  composition,  since  we  do  not  know  the 
height  at  which  the  bulk  of  the  CME  material  originates  (other 
than  that  it  is  coronal). 

These  assumptions  together  result  in  a  lower  limit  for  the 
mass.  The  uncertainty  in  the  CME  mass  becomes  more  signifi¬ 
cant  as  the  central  angle  and/or  spread  of  a  given  CME  departs 
significantly  from  the  plane  of  the  sky.  The  mass  underestima¬ 
tion  is  about  a  factor  of  two  for  CMEs  that  are  <40  from  the 
sky  plane  (Vourlidas  et  al.  2010). 

Other  uncertainties  in  this  procedure  include  exposure  time 
variations  between  event  and  pre-event  images,  improper 
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vignetting  correction,  solar  rotation  effects,  and  the  presence 
of  stars  in  the  field  of  view.  Fortunately,  such  uncertainties  can 
be  minimized  to  a  level  that  is  well  below  that  of  other  fac¬ 
tors  through  proper  calibration  and  careful  choice  of  event  and 
pre-event  images. 

After  obtaining  a  series  of  excess  masses  of  the  CME  as 
a  function  of  time,  we  can  compute  both  the  total  mass  of 
the  CME  and  the  position  and  projected  velocity,  both  for  the 
leading  edge  and  for  the  center  of  mass  of  the  CME.  From  the 
mass  m,  position  /-,  and  velocity  V,  we  can  straightforwardly 
estimate  the  total  kinetic  (C/k  =  (1/2 )mV2)  and  potential 
(Co  =  GMqUi{RqX  —  r — 1 ))  energies  (here  G  is  the  Newtonian 
gravitational  constant  and  MQ  and  Rq  are  the  solar  mass  and 
radius,  respectively).  These  values  are  again  lower  bounds  since 
both  the  mass  and  the  speed  are  projected  quantities.  Vourlidas 
et  al.  (2010)  estimate  that,  for  CMEs  that  are  far  away  from 
the  sky  plane  and  that  have  relatively  small  widths,  the  kinetic 
energy  could  be  as  much  as  a  factor  of  eight  times  larger  than 
the  values  derived  above;  similarly  the  potential  energy  could 
be  as  much  as  twice  as  large  for  such  events.  Elowever,  for  the 
majority  of  events,  the  uncertainties  on  the  quoted  energies  are 
within  a  factor  of  two.  To  obtain  the  kinetic  energy  in  the  solar 
wind  rest  frame  (as  an  estimate  of  the  energy  available  for  shock 
acceleration  of  SEPs;  see  Section  3.6),  we  simply  subtracted 
400  km  s  1  from  the  measured  CME  speed  and  recomputed  the 
kinetic  energy  using  the  speed  in  this  new  reference  frame. 

2.7.  Solar  Energetic  Particles  (SEPs) 

For  the  majority  of  the  events  studied,  it  is  likely  that 
the  interplanetary  SEPs  in  the  events  studied  are  accelerated 
by  CME-driven  shocks.  (A  possible  exception  is  the  2002 
February  20  event,  where  particles  directly  accelerated  in  the 
flare  could  dominate;  see  Chollet  et  al.  2010.)  The  energy 
content  of  the  accelerated  SEPs,  particularly  when  compared 
to  the  kinetic  energy  of  the  CME  in  the  solar  wind  rest  frame, 
is  therefore  an  important  measure  of  the  efficiency  of  SEP 
production  by  the  CME. 

The  energy  content  of  SEPs  that  escape  into  interplanetary 
space  has  been  estimated  by  measuring  the  energy  spectra  of 
electrons  from  ~0.035  to  ~8  MeV,  protons  from  ~0.05  to 
~400  MeV  nucleon-1,  and  abundant  heavier  ions  from  ~0.05 
to  ~  1 00  MeV  nucleon-1.  Estimates  were  made  in  a  number 
of  large  events  from  Solar  Cycle  23  using  a  combination  of 
nine  separate  instruments.  The  proton  spectra  are  based  on 
data  from  the  Ultra-Low  Energy  Isotope  Spectrometer  (ULEIS; 
Mason  et  al.  1998)  and  the  Electron,  Proton,  and  Alpha  Monitor 
(EPAM;  Gold  et  al.  1998)  on  AC£;  from  the  Proton/Electron 
Telescope  (PET;  Cook  et  al.  1993)  on  SAMPEX\  and  from 
the  Energetic  Particle  Sensors  (EPS;  Onsager  et  al.  1996)  on 
NOAA’s  GOES-8  and  GOES-1 1  satellites.  Spectra  of  helium  and 
heavier  ions  were  measured  by  the  Solar  Isotope  Spectrometer 
(SIS;  Stone  et  al.  1998)  on  ACE  and  by  ULEIS.  Also  used 
for  the  2006  events  were  two  STEREO  instruments,  the  Low 
Energy  Telescope  (LET;  Mewaldt  et  al.  2008a)  and  High 
Energy  Telescope  (HET;  von  Rosenvinge  et  al.  2008).  Electron 
measurements  were  provided  by  ACE/ EPAM,  SAMPEX/ PET, 
and  by  the  Electron  Proton  Helium  INstrument  (EPHIN ;  Miiller- 
Mellin  et  al.  1995)  instrument  on  SOHO. 

For  11  of  these  events,  the  energy  spectra  of  H,  He,  and 
abundant  heavier  ions  were  all  fitted  with  common  spectral 
forms  that  include  the  double-power-law  function  of  Band  et  al. 
(1993)  and  the  Ellison  &  Ramaty  (1985)  spectrum — a  power 
law  with  an  exponential  cutoff.  Examples  of  energy  spectra 


and  both  functional  forms  are  given  in  Mewaldt  et  al.  (2005, 
2012)  and  Cohen  et  al.  (2005).  For  the  remainder  of  the  events, 
the  proton  energy  spectra  were  fitted  and  the  contributions  of 
He  and  heavier  ions  were  estimated  using  element  abundances 
measured  for  these  events  by  ULEIS  and  SIS.  The  electron 
contribution  was  measured  in  each  of  the  individual  events  using 
either  EPAM  and  PET,  or  EPAM  and  EPHIN. 

For  all  of  the  fluence  measurements  described  above,  the 
instruments  were  located  near  Earth.  As  in  Emslie  et  al.  (2004), 
we  used  the  measured  near-Earth  fluence  spectra,  typically 
integrated  over  3-5  days,  to  estimate  the  energy  cm-2  that 
escaped  beyond  1  AU  in  the  form  of  SEPs.  To  obtain  this 
estimate,  Emslie  et  al.  (2004)  corrected  for  the  fact  that  SEPs 
can  scatter  back  and  forth  across  1  AU  (providing  multiple 
opportunities  to  be  measured)  using  correction  factors  based  on 
simulations  by  J.  Giacalone  (2002,  private  communication).  A 
similar  approach  was  followed  in  analyzing  the  "Halloween” 
events  (Mewaldt  et  al.  2005)  and  in  a  subsequent  survey  of  17 
events  (Mewaldt  2006).  Mewaldt  et  al.  (2008b)  improved  on 
these  estimates  in  a  study  of  23  SEP  events  from  1997  to  2005 
by  correcting  for  the  fact  that  SEPs  also  lose  energy  as  they 
scatter  on  the  diverging  interplanetary  magnetic  field. 

For  this  work,  we  corrected  for  both  multiple  1  AU  crossings 
and  energy  loss  using  new  simulations  by  Chollet  et  al.  (2010) 
for  four  species  (H,  He,  O,  and  Fe)  with  a  range  of  charge-to- 
mass  ratios.  Chollet  et  al.  (2010)  considered  scattering  mean 
free  paths  X  ranging  from  0.01  to  1  AU,  and  also  varied  the 
radial  and  rigidity  dependence  of  X.  Surprisingly,  the  source 
energy  required  to  account  for  the  accelerated  particles  in  these 
different  scattering  descriptions  varied  by  less  than  a  factor  of 
~2.  This  is  apparently  because  the  scattering  and  energy-loss 
processes  compensate  for  each  other — the  more  particles  scatter 
the  more  often  they  cross  1  AU,  but  they  also  lose  more  energy 
in  the  process.  In  this  paper,  we  have  used  their  form  of  X 
derived  from  quasi-linear  theory  (see  Equation  (3)  in  Chollet 
etal.  2010). 

To  relate  the  measured  near-Earth  values  of  MeV  cm-2 
to  the  integrated  contribution  of  SEPs  escaping  through  a 
1  AU  sphere  surrounding  the  Sun,  we  need  to  know  how 
SEPs  from  a  given  source  location  are  distributed  in  longi¬ 
tude  and  latitude.  Emslie  et  al.  (2004)  assumed  that  the  SEP 
fluence  at  Earth  falls  off  exponentially  with  <?-folding  separa¬ 
tions  of  35°  for  latitude,  45°  for  longitude  in  western  events, 
and  25°  for  longitude  in  eastern  events.  Since  then,  Lario 
et  al.  (2006)  have  measured  the  longitudinal  distribution  of 
SEPs  using  two-  and  three-spacecraft  data  from  the  two  Helios 
spacecraft  and  the  Interplanetary  Monitoring  Platform-8 
(IMP-8).  They  adopted  a  Gaussian  spatial  distribution  given  by 
F  —  Fa  exp[— (O  —  0(|)2/2ct2  |,  where  O  is  the  longitude  of  the 
observer,  0()  is  located  25?  8  east  of  the  point  of  best  solar  wind 
connection  for  a  450  km  s-1  solar  wind  (~W52°),  and  a  =  38°. 
We  use  their  result  for  the  fluence  of  4-13  MeV  protons  and  we 
assume  it  also  applies  to  latitude  differences.  By  using  this  re¬ 
lation  with  the  measured  flare  location  and  the  near-Earth  value 
for  the  escaping  MeV  cm-2,  we  obtained  the  source  energy  re¬ 
quired  to  supply  SEPs  escaping  over  a  1  AU  sphere  centered 
on  the  Sun.  The  results  are  tabulated  in  the  “SEP”  column  of 
Table  1 .  Note  that  we  have  limited  this  study  to  SEP  events  with 
source  regions  ranging  from  E60°  to  W90°  in  longitude.  Be¬ 
yond  this  range  the  Gaussian  longitude  distribution  adopted  by 
Lario  et  al.  (2006)  drops  off  very  rapidly  and  the  longitude  cor¬ 
rections  become  considerably  greater  and  more  uncertain.  The 
typical  uncertainty  in  SEP  energy  is  conservatively  estimated  to 
be  a  factor  of  three. 
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2.8.  Non-potential  Energy  in  the  Active  Region  Magnetic  Field 

It  is  commonly  believed  that  the  fundamental  energy  source 
for  an  SEE  lies  in  current-carrying  magnetic  fields.  In  such 
a  scenario,  the  free  energy  available  to  power  the  event 
is  the  excess  “non-potential”  magnetic  energy — the  energy 
above  the  minimum-energy,  potential  (i.e.,  current-free)  field  to 
which  the  field  can  relax.  The  estimated  available  non-potential 
magnetic  energies  of  the  active  region  producing  the  SEEs  are 
listed  in  the  column  labeled  “Mag”  of  Table  1 .  The  estimates 
were  made  from  full-disk  LOS  magnetograms  obtained  from 
the  Michelson  Doppler  Imager  (MDI;  Scherrer  et  al.  1995)  on 
SOHO,  using  the  method  described  by  Welsch  et  al.  (2009)  and 
outlined  below. 

The  2002  July  20  event  (Event  5)  was  located  behind  the  east 
limb.  Although  its  source  is  therefore  uncertain,  subsequent 
active  region  maps  suggest  that  this  event  probably  occurred  in 
AR  10039,  the  same  source  region  for  the  2002  July  23  event 
(Event  6).  Consequently,  the  estimated  non-potential  magnetic 
energy  for  the  latter  event  is  a  plausible  estimate  for  the  non¬ 
potential  magnetic  energy  in  the  former,  and  indeed  the  use 
of  such  an  estimate  is  consistent  with  the  procedure  used  to 
estimate  the  non-potential  magnetic  energy  content  in  other 
near-limb  events  (e.g.,  W72°  for  2002  February  20  (Event  1) 
and  W84°  for  2002  April  21  (Event  2)),  for  which  magnetic 
field  measurements  obtained  when  the  pertinent  active  region 
was  near  to  disk  center  were  used.  However,  because  there  is 
still  a  finite  possibility  of  misidentification  of  the  active  region, 
we  have  chosen  not  to  use  this  estimated  value  either  in  Table  1 
or  in  the  pertinent  plots  of  Section  3. 

Numerous  efforts  have  been  undertaken  to  estimate  non¬ 
potential  magnetic  energies  in  active  regions  near  disk  center. 
The  methods  include:  (1)  using  the  magnetic  virial  theorem 
estimates  from  chromospheric  vector  magnetograms  (Metcalf 
et  al.  1995,  2005),  (2)  semiempirical  flux-rope  modeling  using 
Ha  and  EUV  images  with  MDI  LOS  magnetograms  (Bobra 
et  al.  2008),  and  (3)  MHD  modeling  (Metcalf  et  al.  1995;  Jiao 
et  al.  1997)  and  non-potential  field  extrapolation  based  upon 
photospheric  vector  magnetograms  (Guo  et  al.  2008;  Schrijver 
et  al.  2008;  Thalmann  &  Wiegelmann  2008;  Thalmann  et  al. 
2008).  These  methods  are  labor  intensive,  and  uncertainties  in 
their  energy  estimates  are  large.  For  example,  error  bars  on  virial 
free-energy  estimates  can  exceed  the  potential  magnetic  energy. 
Also,  there  is  considerable  scatter  in  estimates  from  studies 
that  employ  several  methods  to  analyze  the  same  data  (e.g., 
Schrijver  et  al.  2008).  A  couple  of  generalizations,  however, 
can  be  made.  Free  energies  determined  by  virial  methods 
matched  or  exceeded  the  potential  field  energy,  while  free 
energies  estimated  using  other  techniques  typically  amounted 
to  a  few  tens  of  percent  of  the  potential  field  energy.  Published 
values  for  free  energies  in  analytic  (Schrijver  et  al.  2006)  and 
semiempirical  (Metcalf  et  al.  2008)  fields  meant  to  model  solar 
fields  also  hover  around  a  few  tens  of  percent  of  the  potential 
field  energy. 

We  estimated  the  free  (i.e.,  non-potential)  magnetic  energies 
listed  in  Table  1  to  be  30%  of  the  potential  magnetic  energy 
determined  from  MDI  full-disk  LOS  magnetograms.  This  is 
believed  to  be  a  conservative  estimate  but  it  has  the  advantage 
that  it  can  be  readily  determined  for  most  of  the  events. 
Some  of  the  events  arose  from  limb  active  regions,  for  which 
simultaneous  magnetograms  are  unavailable.  Even  if  vector 
magnetograms  were  available,  uncertainties  in  free  energies 
would  still  be  large.  With  published  virial  free-energy  estimates 
ranging  to  a  few  times  the  potential  energy,  it  is  possible  that 


the  true  free  energy  could  exceed  our  estimates  by  a  factor 
of- 10. 

Apart  from  two  cases  where  flux  was  clearly  emerging 
near  the  time  of  the  event  (Events  29  and  38),  we  calcu¬ 
lated  the  potential  magnetic  energies  from  magnetograms  in 
which  each  event’s  source  active  region  was  near  the  disk’s 
central  meridian,  assuming  a  rigid  rotation  rate  of  13°  day-1. 
This  means  the  energy  estimates  were  sometimes  made  a  few 
days  before  or  after  a  given  event.  Fields  were  assumed  to 
be  radial,  so  each  pixel’s  LOS  field  strength  Blos  was  di¬ 
vided  by  the  cosine  of  the  heliocentric  angle  between  the 
pixel  and  the  sub-observation  point,  to  generate  an  estimated 
radial  field,  Br.  Using  a  Mercator  projection  (Welsch  et  al. 
2009),  the  corrected  pixel  values  were  then  interpolated  onto 
a  two-dimensional  plane.  Next,  the  scalar  potential  /,  where 
B  =  -VX,  was  determined  using  a  Green’s  function  method. 
Finally,  the  magnetic  energy  Um  was  estimated  by  integrating 
('/  Br/Stt)  over  manually  defined  cropping  windows  that  con¬ 
tained  each  active  region.  Images  of  the  magnetograms  used,  as 
well  as  deprojected  data  with  cropping  windows,  are  online 
at  http://solarmu.ri.  ssl.  berkeley.  edu/~welsch/pablic/meetings/ 
SADOSC24/. 

In  several  cases  in  Table  1  (e.g..  Events  22-25),  the  same 
value  of  the  magnetic  energy  is  given  to  adjacent  events  up  to 
five  days  apart  from  the  same  active  region.  This  is  because 
the  magnetic  energy  was  estimated  from  LOS  magnetograms 
taken  when  the  active  region  was  close  to  disk  center.  These 
estimates  become  increasingly  unreliable  as  the  active  region 
moves  away  from  disk  center.  Thus,  although  it  is  very  likely  that 
the  active  region’s  magnetic  fields  evolved  substantially  over  the 
time  between  events,  there  is  no  way  to  reliably  quantify  these 
changes  from  the  available  magnetograms.  This  problem  will 
be  alleviated  with  the  now  regularly  available  vector  magnetic 
field  measurements  from  the  Helioseismic  and  Magnetic  Imager 
(HMI;  Scherrer  et  al.  2012)  on  the  Solar  Dynamics  Observatory 
( SDO ),  which  can  be  used  to  estimate  the  energy  in  regions 
located  away  from  disk  center. 

3.  COMPARISONS  OF  ENERGETIC  COMPONENTS 

In  this  section,  we  present  comparisons  of  the  energy  contents 
of  the  various  components  discussed  above,  through  a  series  of 
figures  each  showing  logarithmic  scatter  plots  (the  “forest”)  of 
the  energy  content  of  one  component  versus  the  energy  content 
of  another,  for  all  events  (“trees”)  for  which  data  are  available 
for  both  selected  components.  The  scatter  of  the  points  around 
the  logarithmic  centroid  is  due  both  to  the  true  range  of  energies 
of  the  selected  events  and  to  the  often  large  uncertainties  (up  to 
2.5  orders  of  magnitude  in  some  cases)  in  the  energy  estimates 
of  each  component.  If  the  uncertainties  are  random,  then  the 
centroid  location  gives  an  indication  of  the  average  ratio  of  the 
energies  of  the  two  components  being  plotted,  and  the  scatter 
of  the  points  about  the  centroid  provides  a  measure  of  the 
overall  uncertainty  in  that  ratio.  Any  “outlier”  point  indicates 
an  anomalous  event,  which  could  simply  identify  an  unusually 
large  or  small  event,  or  which  could  reveal  intrinsic  differences 
in  the  distribution  of  energies  between  the  different  components 
or  some  error  in  the  energy  estimates  for  the  event  in  question. 

The  component  energy  comparisons  are  discussed  in  the  fol¬ 
lowing  subsections,  with  associated  plots  given  in  Figures  1-3. 
All  plots  have  the  same  four  orders-of-magnitude  range  on  each 
axis,  so  that  the  degree  of  spread  in  a  particular  energetic  compo¬ 
nent  can  be  readily  visualized.  In  each  plot,  all  events  that  have 
measured  energies  for  both  components  are  shown.  The  points 
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Thermal  Radiated  Energy  (1030  ergs) 

(a)  GOES  1-8  A  band  vs.  the  total  energy 
radiated  from  the  SXR-emitting  plasma. 


(b)  Total  energy  radiated  from  the  SXR- 
emitting  plasma  vs.  the  peak  thermal 
energy  content  of  that  plasma. 


(c)  Peak  thermal  energy  in  the  SXR- 
emitting  plasma  vs.  energy  in  flare- 
accelerated  nonthermal  particles  (electrons 
plus  ions  when  available). 


(d)  Thermal  radiated  energy  from  the 
SXR-emitting  plasma  vs.  the  energy  in 
flare-accelerated  nonthermal  particles  (elec¬ 
trons  plus  ions  when  available). 


Figure  1.  Scatter  plots  of  different  energy  components,  in  units  of  1030  erg.  Each  plot  includes  all  events  for  which  measurements  are  available  for  both  components. 
The  points  are  indicated  by  triangles,  except  for  the  “outlier”  points  lying  outside  the  2 a  ellipse  (see  below),  which  are  instead  labeled  by  their  event  numbers  (Table  T) 
and  are  located  at  the  center  of  the  respective  numbers.  The  location  of  the  logarithmic  centroid,  defined  by  Equation  (Al),  for  all  the  events  in  the  plot  is  shown  by  a 
bull’s  eye  with  its  X  and  Y  coordinates  listed  in  the  upper  left  corner  of  the  plot.  The  three  diagonal  dashed  lines  represent  the  1%,  10%,  and  100%  ratios  between  the 
plotted  components.  Lines  of  constant  logarithmic  average  event  energy  are  shown  by  tick  marks  every  order  of  magnitude  along  the  lines  of  constant  ratio.  The  major 
and  minor  axes  of  the  ellipse  are  defined  by  ±2  times  the  rms  deviation  of  the  points,  respectively,  parallel  (Equation  (A5))  and  perpendicular  (Equation  (A4))  to  the 
line  of  constant  ratio  passing  through  the  centroid  (see  the  text  and  Appendix  for  discussion).  Points  outside  this  ellipse  are  considered  as  outliers  and  are  discussed  in 
Section  5. 


are  indicated  by  triangles,  except  for  the  “outlier”  points  lying 
outside  the  2a  ellipse  (see  below),  which  are  instead  labeled  by 
their  event  numbers  (Table  1)  and  are  located  at  the  center  of 
the  respective  numbers.  Events  with  only  upper  and/or  lower 
limits  are  generally  not  shown.  However,  we  have  included  per¬ 
tinent  data  for  the  behind-the-limb  Event  5  and  for  Event  25,  for 
which  there  is  some  ambiguity  in  the  CME  kinetic  energy  and 
we  have  hence  used  the  geometric  mean  of  the  two  estimates. 
The  logarithmic  centroid  is  shown  by  a  bull’s  eye  with  its  X 
and  Y  coordinates,  calculated  using  Equation  (Al),  given  in  the 
upper  left  corner  of  the  plot.  The  three  diagonal  dotted  lines 
are  lines  of  constant  ratio  R  as  defined  by  Equation  (A2),  with 
R  —  Y / X  =  100%,  10%,  and  1%.  These  lines  each  have  tick 
marks  showing  the  overall  “size”  of  the  event  A  =  \/XY — see 
Equation  (A3).  The  dashed-line  ellipse  shows  the  ±2er  locus; 
the  widths  of  this  ellipse  perpendicular  (Equation  (A4))  and 
parallel  (Equation  (A5))  to  the  lines  of  constant  ratio  are  mea¬ 
sures  of  the  2 a  spread  in  the  energy  ratio,  R,  and  the  event 
size.  A,  respectively.  Points  outside  this  ellipse  are  considered 
as  “outliers”  and  will  be  discussed  in  Section  5. 


Table  3  lists,  for  each  plot,  the  energetic  components  involved, 
the  value  of  the  logarithmic  centroid  energies  and  their  ratio,  and 
the  root  mean  square  (rms)  spreads  in  the  values  of  the  ratio  R 
and  the  size  A.  Also,  to  quantify  possible  trends  of  one  parameter 
versus  the  other,  Table  3  lists  the  Spearman’s  rank  correlation 
coefficient  p,  a  quantity  that  measures  the  correlation  between 
their  rank  orders  (lowest  -»  highest)  of  the  variables.  The  formal 
equations  used  to  determine  these  different  parameters  are  given 
in  the  Appendix. 

3.1.  Radiated  Energy  in  the  GOES  1-8  A  Band  versus  Radiated 
Energy  from  SXR-emitting  Plasma 

Figure  1  (a)  shows  the  scatter  plot  for  the  radiated  energy  in  the 
GOES  1-8  A  band  versus  the  total  energy  radiated  from  the  hot 
SXR-emitting  plasma.  The  points  are  closely  rank-correlated 
(p  =  0.96)  and  also  cluster  very  closely  in  the  perpendicular 
(R)  direction,  showing  that  the  energy  radiated  in  the  GOES 
1-8  A  band  is  a  relatively  constant  fraction  ( R  =  0.05)  of  the 
total  energy  radiated  from  the  SXR-emitting  plasma.  Indeed, 
a  regression  analysis  of  the  data  in  Table  1  shows  that  the 


9 


The  Astrophysical  Journal,  759:71  (18pp),  2012  November  1 


Emslie  et  al. 


(a)  Energy  in  flare-accelerated  ions  vs.  en¬ 
ergy  in  flare-accelerated  electrons. 


(b)  Energy  in  SEPs  vs.  CME  kinetic  energy 
in  the  rest  frame  of  the  solar  wind. 


CO 


(c)  Energy  in  SEPs  vs.  energy  in  flare-  (d)  Bolometric  radiated  energy  vs.  non¬ 
accelerated  ions.  potential  magnetic  energy  in  the  active  region. 

Figure  2.  Same  as  Figure  1  for  different  combinations  of  energy  components,  as  indicated  on  the  axis  labels. 


Table  3 

Parameters  Represented  in  Scatter  Plots 


Figure 

No. 

Plotted  Components 

Log.  Centroida 

Rb 

rmsc 

Pd 

No.  of 

Events 

X-Axis 

7-Axis 

X 

(1030  erg) 

Y 

(1030  erg) 

R 

A 

1(a) 

Rad.  from  hot  plasma 

GOES  1-8  A 

12 

0.6 

0.05 

0.17 

0.51 

0.96 

38 

1(b) 

Peak  thermal  energy 

Rad.  from  hot  plasma 

3.9 

11 

2.7 

0.22 

0.34 

0.82 

26 

1(c) 

Electrons+ions 

Peak  thermal  energy 

34 

3.9 

0.11 

0.43 

0.31 

0.36 

26 

1(d) 

Electrons+ions 

Rad.  from  hot  plasma 

34 

11 

0.31 

0.43 

0.34 

0.46 

26 

2(a) 

Electrons 

Ions 

32 

11 

0.34 

0.63 

0.52 

0.45 

14 

2(b) 

CME  KE  (SW  frame) 

SEP 

110 

4.0 

0.04 

0.49 

0.47 

0.47 

20 

2(c) 

Ions 

SEP 

23 

6.2 

0.27 

0.38 

0.35 

0.20 

6 

2(d) 

Magnetic 

Bolometric 

890 

55 

0.06 

0.43 

0.38 

0.56 

37 

3(a) 

Magnetic 

CME  KE+PE 

1000 

200 

0.19 

0.39 

0.34 

0.68 

23 

3(b) 

CME  KE+PE 

Bolometric 

200 

71 

0.35 

0.43 

0.38 

0.54 

24 

3(c) 

Bolometric 

Electrons+ions 

49 

34 

0.71 

0.47 

0.35 

0.37 

26 

3(d) 

Bolometric 

Rad.  from  hot  plasma 

57 

12 

0.21 

0.21 

0.47 

0.92 

38 

Notes. 

a  X  and  Y  values  of  the  logarithmic  centroid,  computed  using  Equation  (Al). 
b  R  =  Y /X,  the  ratio  of  Y  and  X  values  of  the  logarithmic  centroid  computed  using  Equation  (A2). 

c  rms  (root  mean  square)  values  of  R  =  Y/X  and  A  =  y/XY,  computed  using  Equations  (A4)  and  (A5).  The  rms  values  of  R  and  A,  respectively,  quantify  the  scatter 
perpendicular  and  parallel  to  the  line  of  constant  energy  ratio  that  passes  through  the  logarithmic  centroid. 

d  Spearman’s  rank  correlation  coefficient — a  non-parametric  measure  of  statistical  dependence  between  two  variables — see  Equation  (A6). 
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cn 

QJ 


(a)  CME  total  energy  (kinetic  +  poten¬ 
tial,  in  the  rest  frame  of  the  Sun)  vs.  the 
free  (non-potential)  magnetic  energy  of  the 
active  region. 


Bolometric  Rodioted  Energy  (1030  ergs) 

(c)  Energy  in  flare-accelerated  nonthermal 
particles  (electrons  and  ions)  vs.  bolometric 
radiated  energy. 

Figure  3.  Same  as  Figure  1  for  different  combinations 

ratio  (best  estimate  ±ler)  of  the  total  energy  radiated  by 
SXR-emitting  thermal  plasma  in  the  flare  to  the  observed  GOES 
1-8  A  flux  is  15. 4 ±0.8.  This  strong  correlation  is  not  surprising 
since  both  energy  components  plotted  are  calculated  from  the 
GOES  SXR  fluxes;  the  scatter  about  the  trend  line  arises  from 
the  differences  in  the  temperatures  of  the  different  events. 

3.2.  Thermal  Radiated  Energy  versus  Peak  Thermal  Energy 

Figure  1  (b)  shows  the  scatter  plot  of  the  total  energy  radiated 
from  hot  SXR-emitting  plasma  versus  the  peak  thermal  energy 
content  of  that  plasma.  The  relatively  tight  correlation  ( p  — 
0.82)  between  these  two  components  is  expected,  since  both 
parameters  refer  to  the  same  SXR-emitting  plasma.  Event  30 
(the  M6.4  event  on  2005  August  25)  is  the  most  extreme  outlier 
but  it  is  almost  equally  weak  in  both  energy  components.  On 
average,  the  total  energy  radiated  exceeds  the  peak  thermal 
energy  content  by  a  factor  of  ~3  (R  =  2.7  in  Table  3),  implying 
continuous  re-energization  of  the  SXR-emitting  material  as  the 
flare  progresses. 

3.3.  Peak  Thermal  Energy  versus  Energy  in  Flare-accelerated 
Nonthermal  Particles 

Figure  1(c)  shows  the  scatter  plot  for  the  peak  thermal  energy 
in  the  hot  SXR-emitting  plasma  versus  the  energy  in  flare- 
accelerated  nonthermal  particles  (electrons  plus  ions,  when 


(b)  Bolometric  radiated  energy  vs.  the 
CME  total  energy  (kinetic  +  potential,  in 
the  rest  frame  of  the  Sun). 


(d)  Total  energy  radiated  by  SXR-emitting 
plasma  vs.  bolometric  radiated  energy. 


of  energy  components,  as  indicated  on  the  axis  labels. 

available).  There  is  substantially  greater  spread  in  the  points 
compared  to  Figure  1(b),  but  nevertheless  a  reasonable  bunching 
of  the  points.  The  maximum  spread  of  the  points  is  less  than  two 
orders  of  magnitude  in  either  parameter.  On  average,  the  energy 
in  the  flare-accelerated  nonthermal  particles  exceeds  the  peak 
thermal  energy  by  almost  an  order  of  magnitude  (R  =  0.11), 
indicating  that  there  is  easily  sufficient  power  in  the  particles 
to  create  the  SXR-emitting  thermal  plasma.  This  conclusion 
is  reinforced  by  the  fact  that  the  energy  in  flare-accelerated 
electrons  is  a  lower  limit  (see  Section  2.4)  and  is  in  agreement 
with  earlier  comparisons  of  flare-accelerated  electrons  versus 
peak  thermal  energy — see,  e.g.,  de  Jager  et  al.  (1986)  and  Saint- 
Hilaire  &  Benz  (2002). 

3.4.  Thermal  Radiated  Energy  versus  Energy  in 
Flare-accelerated  Nonthermal  Particles 

Figure  1  (d)  shows  the  scatter  plot  for  the  total  energy  radiated 
by  the  SXR-emitting  thermal  plasma  versus  the  energy  in 
flare-accelerated  nonthermal  particles  (electrons  and  ions  when 
available).  This  figure  combines  information  already  evident  in 
Figures  1(b)  and  (c).  It  shows  that  the  energy  in  accelerated 
electrons  and  ions  during  a  flare  is  not  only  sufficient  to  supply 
the  peak  energy  of  the  SXR-emitting  plasma  (Figure  1(c)), 
but  it  is  also  high  enough  (by  a  factor  of  ~3,  R  —  0.31) 
to  account  for  the  radiation  from  this  plasma  throughout  the 
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event.  As  discussed  in  Section  2,  it  follows  that  a  significant 
fraction  of  the  energy  in  flare-accelerated  nonthermal  particles  is 
deposited  by  thermal  conduction  into  lower-temperature  plasma 
and  ultimately  radiated  in  optical  and  EUV  wavebands  (see 
Emslie  et  al.  2005).  Again,  Event  30  (the  M6.4  event  on  2005 
August  25)  is  the  only  “outlier”  in  this  plot,  reflecting  the  low 
values  of  both  the  thermal  and  nonthermal  energy  components. 

3.5.  Flare-accelerated  Ions  versus  Electrons 

Figure  2(a)  shows  the  scatter  plot  for  the  energy  in  flare- 
accelerated  ions,  as  determined  from  the  RHESSI  gamma-ray 
observations  (Section  2.5),  versus  the  energy  in  flare-accelerated 
electrons,  as  determined  from  RHESSI  hard  X-ray  observations 
(Section  2.4). 

As  pointed  out  in  Section  2.4,  the  energy  in  electrons  is 
critically  dependent  on  the  low-energy  cutoff,  Enm],  that  is 
assumed  for  the  electron  spectrum.  Since  the  largest  value  of 
Zimin  that  gives  an  acceptable  fit  to  the  data  is  used  for  each 
spectrum,  the  total  electron  energy  values  are  lower  limits  with 
order-of-magnitude  uncertainties.  As  explained  in  Section  2.5, 
the  situation  for  the  ion  energies  is  even  worse,  both  because 
of  the  spread  in  the  observed  2.223  MeV  fluences  and  because 
of  the  need  to  extrapolate  the  ion  flux  at  energies  above  30  MeV, 
as  derived  from  these  2.223  MeV  line  fluences,  to  the  ion  flux 
above  1  MeV.  Because  of  these  large  uncertainties  in  both  the 
electron  and  ion  energies,  there  is  a  much  wider  scatter  than  in 
the  plots  in  Figure  1.  However,  with  some  notable  exceptions 
that  have  almost  two  orders  of  magnitude  more  energy  in  the 
electrons  than  in  the  ions  (Events  9  and  15),  the  electron  and  ion 
energies  are  generally  comparable  within  an  order  of  magnitude. 
This  result  is  in  agreement  with  the  claims  by  Ramaty  et  al. 
(1995)  and  Ramaty  &  Mandzhavidze  (2000)  and  has  significant 
consequences  for  particle  acceleration  models. 

3.6.  SEP  Energy  versus  CME  Kinetic  Energy 

in  the  Solar  Wind  Rest  Frame 

Figure  2(b)  shows  the  scatter  plot  for  the  energy  in  the 
accelerated  SEP  population  versus  the  kinetic  energy  of  the 
CME  in  the  rest  frame  of  the  solar  wind.  We  use  the  solar  wind 
rest  frame  since  a  shock  can  be  formed  and  SEPs  accelerated 
only  if  the  CME  is  traveling  at  least  as  fast  as  the  solar  wind 
speed  (Mewaldt  et  al.  2008a).  Lacking  knowledge  of  the  solar 
wind  speed  low  in  the  corona  for  each  event,  we  have  simply 
subtracted  400  km  s' 1  from  the  measured  CME  speed  in  order 
to  estimate  the  kinetic  energy  available  for  accelerating  particles 
via  shock  acceleration. 

Most  of  the  SEP  values  cluster  between  1%  and  10%  of 
the  CME  kinetic  energy.  Comparing  the  nine  events  that  are 
common  to  both  Mewaldt  et  al.  (2008b)  and  this  study  (Events 
2,  4,  7,  8,  11,  12,  13,  14,  and  16),  the  SEP/CME  ratio  was 
5.8%  in  Mewaldt  et  al.  (2008b)  and  is  4%  here  (R  —  0.04  in 
Table  3).  Overall,  as  a  result  of  several  changes  in  the  analysis, 
the  SEP  energy  estimates  in  this  paper  are  reduced  from  those 
in  Mewaldt  et  al.  (2005,  2008b)  by  an  average  of  40%.  One  of 
these  changes  is  the  adoption  of  the  longitude  correction  of  Lario 
et  al.  (2006),  which  results  in  changes  of  as  much  as  a  factor 
of  two  in  the  energy  estimates  of  individual  SEP  events.  For 
the  nine  events  in  common  with  this  paper  and  Mewaldt  et  al. 
(2008b),  the  differences  due  to  longitude/latitude  corrections 
alone  ranged  from  —52%  to  +51%  with  a  mean  difference  of 
— 16%;  for  all  20  events  shown  in  Figure  2(b),  the  average  effect 
is  a  ~10%  decrease  in  energy.  Another  change  is  the  adoption 


of  new  corrections  for  particles  crossing  1  AU  multiple  times 
and  for  adiabatic  energy  loss  that  are  based  on  the  energy-  and 
species-dependent  simulations  of  Chollet  et  al.  (2010);  these 
changes  resulted  in  a  further  decrease  averaging  ~30%  in  the 
SEP  energy  content  of  the  nine  events  in  common  with  Mewaldt 
et  al.  (2008b).  Finally,  the  SEP  spectra  in  this  study  were 
integrated  from  0.03  MeV  nucleon'1  to  300  MeV  nucleon-1, 
rather  than  from  0.01  MeV  nucleon-1  to  1000  MeV  nucleon-1 
as  in  Emslie  et  al.  (2004)  and  Mewaldt  et  al.  (2005,  2008b). 
The  increase  in  the  low-energy  limit  to  0.03  MeV  nucleon-1 
represents  a  more  realistic  threshold  for  injection  into  the  shock 
acceleration  process  (see,  e.g.,  Li  et  al.  2012),  and  results  in  a 
typical  reduction  in  the  energy  content  by  ~5%.  The  change  in 
the  upper  limit  has  a  negligible  effect. 

Overall,  this  new  analysis  confirms  that  the  SEP  energy  is  a 
small,  but  not  insignificant,  fraction  of  the  CME  kinetic  energy 
in  most  large  events. 

3. 7.  SEP  Energy  versus  Energy  in  Flare-accelerated  Ions 

Figure  2(c)  shows  the  scatter  plot  for  the  total  energy  in 
SEPs  versus  the  energy  in  flare-accelerated  ions,  as  determined 
from  RHESSI  gamma-ray  observations.  There  are  only  a  limited 
number  of  events  that  can  be  compared,  but  in  those  few  cases 
there  is  comparable  energy  in  the  ions  and  SEPs  (R  =  0.27). 
At  first  sight,  this  appears  to  conflict  with  the  study  by  R.  A. 
Mewaldt  (2012,  in  preparation)  that  shows  the  number  of 
>30  MeV  SEP  protons  measured  in  interplanetary  space  is 
generally  much  higher  than  the  number  of  >30  MeV  protons 
interacting  in  the  solar  atmosphere  (Shih  et  al.  2009).  However, 
this  difference  can  be  accounted  for  in  a  number  of  ways.  First, 
as  noted  in  Section  2.5,  the  inferred  energy  in  flare-accelerated 
ions  is  based  on  a  very  uncertain  spectral  extrapolation  over 
more  than  an  order  of  magnitude  in  proton  energy  (from  30  MeV 
down  to  1  MeV).  To  obtain  the  energy  in  flare-accelerated  ions 
(Section  2.5),  we  have  assumed  a  spectral  index  of  4,  which 
results  in  a  significantly  higher  total  energy  than  would  be 
obtained  if  the  extrapolation  was  performed  with  the  much  lower 
spectral  indices  representative  of  SEP  spectra  measured  in  situ 
at  1  AU.  Further,  for  the  well-observed  2003  October  28  flare 
(Event  12),  we  obtained  an  ion  spectral  index  of  3.4  at  energies 
between  ~3  and  50  MeV,  so  that  the  energy  content  in  ions  could 
be  significantly  lower  than  we  have  used  here,  especially  if  the 
spectrum  hardens  even  more  between,  say,  1  and  20  MeV,  a 
feature  that  is  clearly  seen  in  SEP  spectra.  Second,  SEP  protons 
typically  carry  >80%  of  the  SEP  ion  energy,  whereas  flare- 
accelerated  protons  carry  only  about  one-third  of  the  ion  energy. 

3.8.  Bolometric  Radiated  Energy  versus  Magnetic  Energy 

Figure  2(d)  shows  the  scatter  plot  for  the  bolometric  radiated 
energy  versus  the  non-potential  (free)  magnetic  energy  in  the 
active  region.  Here,  we  have  a  relatively  tight  bunching  with 
little  more  than  1 .5  orders  of  magnitude  range  in  each  parameter. 
It  should  be  emphasized  that  the  plotted  bolometric  radiated 
energies  are,  with  the  exception  of  the  five  events  noted  in 
Table  1,  not  directly  measured  but  rather  estimates  made  using 
the  FISM  model  (Chamberlin  et  al.  2007,  2008),  and  the 
magnetic  energy  is  only  good  to  an  order  of  magnitude.  With 
this  proviso,  we  find  that  the  average  bolometric  radiated  energy 
is  ~6%  of  the  free  magnetic  energy  ( R  =  0.06),  and  in  all  cases 
the  available  magnetic  energy  exceeds  the  bolometric  radiated 
energy  by  at  least  half  an  order  of  magnitude.  This  is  consistent 
with  the  well-accepted  notion  that  the  reservoir  of  magnetic 
energy  is  sufficient  to  power  the  main  components  of  the  flare. 
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3. 9.  CME  Energy  versus  Magnetic  Energy 

Figure  3(a)  shows  the  scatter  plot  for  the  CME  total  energy 
(potential  +  kinetic)  in  the  rest  frame  of  the  Sun  versus  the 
non-potential  energy  in  the  magnetic  field.  The  CME  energy  is, 
on  average,  only  a  small  fraction  ( R  —  0.19)  of  the  magnetic 
energy,  similar  to  the  value  R  =  0.06  found  for  the  ratio  of 
bolometric  radiated  energy  to  magnetic  energy  (Section  3.8; 
Figure  2(d)).  While  bearing  in  mind  the  very  approximate  values 
of  the  latter,  it  nevertheless  appears,  from  the  results  of  this  and 
the  previous  subsection,  that  much  of  the  available  magnetic 
energy  (some  two-thirds)  is  retained  in  the  active  region  (i.e., 
the  field  does  not  return  to  a  fully  potential  state),  even  after 
the  flare  and  the  ejection  of  the  CME.  This  result  is  consistent 
with  the  generous  limits  established  by  Moore  et  al.  (2012)  on 
the  possible  free  energy  that  an  active  region  magnetic  field  can 
hold  before  it  erupts. 

3.10.  Bolometric  Radiated  Energy  versus  CME  Energy 

Figures  2(d)  and  3(a)  show,  respectively,  that  the  available 
magnetic  energy  is  about  15  times  the  bolometrically  radiated 
energy  in  the  flare  (R  —  0.06)  and  about  5  times  the  CME  energy 
(R  —  0.19).  Given  the  substantial  overlap  of  events  common 
to  both  plots  (Tables  1  and  3),  this  indicates  that,  on  average, 
the  energy  in  the  CME  is  larger  than  the  energy  radiated  by  a 
factor  of  ~3.  Figure  3(b)  confirms  this  result  by  showing  the 
scatter  plot  for  bolometric  radiated  energy  versus  the  CME  total 
energy  (kinetic  +  potential).  On  average,  the  bolometric  energy 
is  indeed  about  half  an  order  of  magnitude  less  than  the  CME 
energy  (R  =  0.35). 

3.11.  Flare-accelerated  Particle  Energy  versus 
Bolometric  Radiated  Energy 

Figure  3(c)  shows  the  scatter  plot  for  the  total  energy  in  flare- 
accelerated  particles  (electrons  plus  ions)  versus  the  bolometric 
radiated  energy.  This  figure  shows  that  the  energy  in  accelerated 
particles  during  a  flare  is  comparable  to  the  total  bolometric 
radiated  energy  from  the  flare  (R  —  0.71),  with  the  ratio  being 
greater  than  unity  in  some  events  and  less  than  unity  in  others.  It 
must  be  recalled  that  while  the  bolometric  radiated  energies 
are  accurate  to  within  a  factor  of  ~2  to  3,  the  energies  in 
flare-accelerated  particles  are  uncertain  to  at  least  an  order  of 
magnitude.  The  energies  in  electrons  are  most  probably  lower 
limits  and  may  well  underestimate  the  true  energy  content  by  up 
to  an  order  of  magnitude.  The  energies  in  ions  may,  however,  be 
overestimates,  depending  on  the  power-law  index  used  for  the 
spectral  extrapolation.  We  can  tentatively  conclude,  however, 
that  there  is  sufficient  energy  in  the  flare-accelerated  particles 
to  account  for  all  the  energy  radiated  in  the  flare.  However, 
this  conclusion  must  somehow  be  verified  by  more  accurate 
estimates  of  the  energy  in  the  flare-accelerated  electrons  and 
ions. 

3.12.  Energy  Radiated  by  Soft  X-Ray  Emitting  Plasma 
versus  Bolometric  Energy 

Figure  3(d)  shows  the  scatter  plot  for  the  total  energy  radiated 
by  SXR-emitting  plasma  versus  the  bolometric  radiated  energy. 
The  relatively  tight  correlation  ( p  =  0.92)  is  to  a  large  extent 
due  to  the  fact  that  most  of  the  bolometric  radiant  energies 
were  computed  using  the  FISM  model,  which  uses  the  radiated 
energy  from  SXR-emitting  plasma  to  estimate  the  bolometric 
energy.  Nevertheless,  the  data  show  that  only  about  one-fifth 
of  the  bolometric  energy  radiated  by  a  flare  is  radiated  by  the 
SXR-emitting  plasma  (R  —  0.21). 
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4.  DISCUSSION 

The  comparisons  of  the  energetics  of  the  different  compo¬ 
nents  shown  in  the  scatter  plots  of  Figures  1-3  are  summarized 
in  Table  3. 

In  attempting  to  draw  any  definitive  conclusions  from  these 
comparisons,  we  must  keep  in  mind  the  limitations  of  our 
analysis.  The  events  in  our  list  cover  the  period  from  RHESSI’s 
launch  in  February  2002  through  2006.  They  include  9  of  the 
1 1  X5  or  greater  events  that  occurred  in  this  time  period.  They 
also  include  the  largest  SEP  events  and  all  significant  RHESSI 
gamma-ray  line  events,  and  hence  represent  events  where  a 
significant  number  of  ions  were  accelerated  to  high  energies, 
either  in  the  flare,  at  the  CME  shock,  or  at  both  locations. 
Most  of  the  CMEs  that  we  have  included  have  kinetic  energies 
of  >  I  (f2  erg.  According  to  Figure  8  of  Gopalswamy  et  al. 
(2004),  only  16  CMEs  observed  from  1997  to  2002  had  kinetic 
energy  >1.4  x  1032  erg.  Adding  in  the  five  CMEs  with  kinetic 
energy  >1.4  x  1032  erg  that  occurred  in  2003  (Events  11,  12, 
13,  14,  and  16  in  Table  1),  we  see  that  only  21  CMEs  with 
kinetic  energy  >1032  erg  occurred  during  1997-2003.  For  this 
period,  Gopalswamy  (2006)  reports  a  total  of  4133  CMEs,  with 
average  kinetic  energy  5  x  1029  erg.  Therefore,  if  we  consider 
the  1997-2003  period  as  typical,  only  21/4133  ~  0.5%  of  all 
CMEs  have  kinetic  energies  as  large  as  those  considered  in  this 
paper.  In  summary,  the  events  studied  represent  the  largest  SEEs 
that  occurred  during  the  period  in  question. 

Despite  the  relatively  large  number  of  events  in  our  list, 
the  range  of  energies  is  typically  only  about  two  orders  of 
magnitude  in  any  of  the  energy  components.  The  uncertainties 
on  the  energy  estimates  of  each  component  are  generally  at 
least  an  order  of  magnitude,  except  for  the  few  cases  where  the 
bolometric  radiated  energy  is  measured  with  an  accuracy  of  a 
factor  of  ~2.  Thus,  in  most  cases,  we  cannot  expect  to  see  any 
significant  trends  with  the  size  of  the  events  over  the  limited 
range  of  our  selected  events.  In  spite  of  these  selection  effects 
and  measurement  limitations,  the  scatter  plots  nevertheless 
reveal  several  useful  results  regarding  large  SEEs  in  general 
and  about  specific  events  in  particular.  With  a  few  notable 
exceptions,  the  points  in  each  plot  are  bunched  together  within 
the  expected  order-of-magnitude  uncertainties.  However,  there 
is  substantially  more  scatter  of  the  points  in  some  parameters 
than  others.  Only  a  few  events  stand  out  as  outliers  in  certain 
plots;  these  outliers  are  discussed  in  Section  5. 

The  general  bunching  of  the  data  points  in  each  scatter  plot 
is  characterized  by  the  logarithmic  rms  deviations  parallel  and 
perpendicular  to  the  line  of  constant  ratio  that  passes  through 
the  logarithmic  centroid.  The  X  and  Y  values  of  the  centroids 
are  shown  on  each  plot  and  listed  along  with  the  rms  values  in 
Table  3.  The  fact  that  the  data  points  generally  bunch  together 
can  be  interpreted  as  an  extension  to  SEEs  of  the  “ big  flare 
syndrome”  (BFS),  a  phrase  coined  by  Kahler  (1982)  based  on 
the  strong  correlations  between  proton  fluxes  and  associated 
microwave  and  hard  X-ray  burst  parameters,  and  a  concept 
which  has  since  come  to  mean  that  each  flare  component  scales 
roughly  linearly  with  some  absolute  measure  of  flare  “size.” 
At  that  time,  before  the  so-called  solar  flare  myth  was  exposed 
(Gosling  1993),  it  was  not  clear  that  SEPs  were  generally  more 
likely  to  be  accelerated  at  CME  shock  fronts,  but  now  the 
same  concept  can  be  applied  to  include  CMEs  with  flares  and 
SEPs.  From  our  data,  it  is  clear  that,  with  some  caveats,  the 
BFS  concept  can  be  applied  to  all  energetic  SEE  phenomena, 
including  the  CME  energy. 
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Notable  exceptions  to  the  scatter  of  the  points  being  consistent 
with  the  expected  uncertainties  in  each  parameter  are  the 
following. 

1.  The  plot  of  the  GOES  1-8  A  integrated  energy  versus  the 
total  energy  radiated  from  the  SXR-emitting  plasma  over 
all  wavelengths  (Figure  1(a))  shows  that  these  two  param¬ 
eters  are  well  correlated  (the  Spearman’s  rank  correlation 
coefficient  p  given  in  Table  3  is  0.96,  the  highest  for  any 
pair  of  parameters).  This  is  not  surprising  since  the  measure¬ 
ments  of  the  two  parameters  are  not  independent — they  both 
use  the  same  GOES  X-ray  data.  The  scatter  of  the  points 
thus  reflects  only  the  range  of  flare  sizes  and  the  different 
temperatures — the  scatter  perpendicular  to  the  line  of  con¬ 
stant  ratio,  only  about  half  an  order  of  magnitude,  is  the 
result  of  different  temperatures,  while  the  scatter  parallel 
to  that  line  is  almost  two  orders  of  magnitude,  reflect¬ 
ing  the  range  of  flare  intensities.  A  regression  analysis  of 
the  data  in  Table  1  shows  that  the  ratio  (best  estimate  ± 
standard  error)  of  the  total  energy  radiated  by  SXR-emitting 
thermal  plasma  in  the  flare  to  the  observed  GOES  1-8  A 
flux  is  15.4  ±  0.8,  i.e.,  that  the  energy  radiated  in  the 
GOES  1-8  A  band  is  about  one-fifteenth  to  one-twentieth 
( R  =  0.05)  of  the  total  energy  radiated  by  the  SXR-emitting 
plasma  over  all  wavelengths.  The  presence  of  the  outlier 
points  for  Events  3  and  8,  several  rms  values  away  from 
the  line  of  constant  ratio,  suggests  that  these  two  events 
are  different  in  that  their  temperature  is  lower  than  the  av¬ 
erage.  This  is  consistent  with  the  conclusion  reached  by 
Feldman  et  al.  (1996)  and  Garcia  (2004)  that  lower  temper¬ 
atures  are  generally  associated  with  smaller  X-ray  peaks. 

2.  The  range  of  the  ion  energies — about  three  orders  of 
magnitude — is  larger  than  for  all  other  parameters.  This 
range  is  especially  evident  in  Figure  2(a),  which  shows 
the  scatter  plot  for  the  energy  in  flare-accelerated  ions 
as  determined  from  the  RIIESSI  gamma-ray  observa¬ 
tions  (Section  2.5)  versus  the  energy  in  flare-accelerated 
electrons  as  determined  from  RHESSI  hard  X-ray  ob¬ 
servations  (Section  2.4).  Since  the  ion  energy  content 
above  1  MeV  was  deduced  by  applying  the  same  8  — 
4  spectral  extrapolation  to  the  energy  content  above 
30  MeV  in  all  events,  this  large  range  in  energy  con¬ 
tents  above  1  MeV  mirrors  exactly  the  spread  in  the 
energy  contents  above  30  MeV  and  hence  the  spread 
in  the  observed  2.223  MeV  line  fluences.  However,  it  must 
be  noted  that  uncertainty  in  the  value  of  <5  results  in  a  further 
large  uncertainty  in  the  ion  energy  content  above  1  MeV 
for  any  specific  event,  since  this  quantity  is  derived  through 
spectral  extrapolation  over  one-and-a-half  orders  of  magni¬ 
tude  in  ion  energy  (Section  2.5).  This  uncertainty  would  act 
to  increase  the  scatter  in  the  energy  content  above  1  MeV 
if  the  value  of  S  was  positively  correlated  with  the  value 
of  the  >30  MeV  energy  content;  alternatively,  it  would 
act  to  decrease  the  scatter  in  energy  content  >1  MeV  if 
the  value  of  <5  was  negatively  correlated  with  the  value  of 
the  >30  MeV  energy  content.  In  this  context,  it  should  be 
noted  that  Shih  et  al.  (2009)  found  a  strong  correlation,  over 
more  than  three  orders  of  magnitude,  between  the  energy 
in  flare-accelerated  ions  and  the  energy  in  electrons  above 
300  keV  (rather  than  the  ~20  keV  lower  cutoff  energy 
used  here).  This  suggests  that  steep  (shallow)  spectra  (high 
(low)  values  of  8)  are  associated  with  low  (high)  values 
of  the  energy  content  at  high  energies,  and  therefore  that 


use  of  individual  spectral  indices  to  create  more  accurate 
energy  estimates  of  the  ion  energy  content  above  1  MeV 
might  reduce  the  scatter  in  the  plot. 

5.  DISCUSSION  OF  OUTLIER  DATA  POINTS 

While  most  of  the  events  in  Table  1  lie  (by  definition)  within 
the  2cr  ellipses  in  the  various  cross-correlation  plots,  there  are 
a  few  notable  exceptions.  We  now  discuss  these  “outliers”  and 
the  possible  reasons  for  their  unusual  energetic  partitioning. 

1.  Event  1.  This  M5.1  event,  on  2002  February  20,  is  one  of 
the  two  events  with  a  relatively  low  ratio  of  CME  energy 
(kinetic  +  potential)  to  magnetic  energy  (Figure  3(a)). 
Table  1  shows  a  paucity  of  data  for  other  energetic  compo¬ 
nents.  Data  for  GOES  1-8  A  emission  and  total  radiated  en¬ 
ergy  from  the  SXR-emitting  plasma  are  available;  the  event 
does  fall  outside  the  2a  ellipse  in  Figure  1(a),  but  only  by 
virtue  of  its  overall  weakness,  not  the  ratio  of  GOES  1-8  A 
emission  to  total  SXR-emitting  energy.  We  believe  that  this 
is  therefore  simply  a  weak  event,  in  which  only  a  small 
fraction  of  the  available  magnetic  energy  was  dissipated. 

2.  Event  3.  This  C5.0  event,  on  2002  May  22,  has  a  rela¬ 
tively  low  (<1%)  ratio  of  GOES  1-8  A  emission  to  total 
SXR-emitting  energy  (Figure  1(a)).  An  event  with  limited 
overall  information  (Table  1),  it  does,  however,  appear  in 
Figure  2(b),  where  a  normal  ratio  of  SEP  to  CME  energy 
is  evident.  It  should  be  noted  that  this  is  the  only  GOES 
C-class  event  in  Table  1  and  we  therefore  simply  categorize 
this  event  as  a  weak  GOES  event,  possibly  due  to  the  low 
temperature  (<9  MK)  of  the  SXR-emitting  plasma.  This 
event  is  the  third  in  a  sequence  of  events  starting  with  an 
Ml. 5  event  peaking  at  21:29  UT  on  May  21,  followed  by 
a  C9.7  event  at  00:30  on  May  22,  and  the  C5.0  event  in 
question  at  03:34  UT.  RHESSI  saw  parts  of  each  of  these 
events  and  the  6—12  keV  images  show  that  they  came  from 
three  distinctly  different  locations  with  the  following  spa¬ 
tial  centroid  coordinates  (inarcseconds):  (—550, 270),  (880, 
—330),  and  (750,  —350),  respectively.  The  peak  tempera¬ 
tures  derived  from  the  GOES  data  for  the  three  events  are 
13,  11,  and  8  MK,  respectively.  Thus,  this  event  was  much 
cooler  than  the  other  larger  events  in  our  list. 

3.  Event  8.  Data  for  this  M4.6  event,  on  2002  November  9, 
are  available  for  all  energetic  components  other  than  flare- 
accelerated  ions  (Table  1);  consequently,  the  event  appears 
in  most  of  the  scatter  plots  in  Figures  1-3.  In  most  of  the 
plots,  the  event  is  situated  within  the  general  bunching  of 
points.  However,  Figure  3(b)  reveals  a  low  ratio  of  bolo- 
metric  energy  to  CME  energy,  while  Figure  3(c)  shows  a 
similarly  high  ratio  of  energy  in  accelerated  particles  to 
bolometric  radiated  energy.  Together,  these  point  simply 
to  an  event  with  a  relatively  low  bolometric  radiance,  as 
inferred  from  the  FISM  model  (Section  2.2).  This  arises 
because  of  the  relatively  short  duration  of  the  GOES  event 
(13  minutes  from  peak  1-8  A  flux  to  50%  of  peak).  As 
pointed  out  by  Chamberlin  et  al.  (2012),  “the  total  radi¬ 
ated  output  of  flares  depends  more  on  the  flare  duration 
than  the  typical  GOES  X-ray  peak  magnitude  classifica¬ 
tion.”  The  relatively  low  temperature  of  this  flare  (peak 
value  of  13  MK)  also  explains  the  relatively  low  ratio  of 
GOES  1-8  A  integrated  flux  to  thermal  radiated  energy 
(Figure  1(a)).  This  event  also  appears  in  Figure  2(b) 
as  having  a  marginally  low  ratio  of  SEP  energy  to 
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CME  kinetic  energy,  but  this  is  a  presumably  unrelated 
phenomenon. 

4.  Event  9.  This  X1.4  event,  on  2003  May  27,  has  the  low¬ 
est  ratio  of  flare-accelerated  ion  energy  to  flare-accelerated 
electron  energy  (Figure  2(a)).  The  ratios  of  peak  thermal  en¬ 
ergy  and  broadband  SXR  radiated  energy  to  energy  in  flare- 
accelerated  nonthermal  particles  are  average  (Figures  1(c) 
and  (d)),  and  the  ratio  of  GOES  1-8  A  to  total  SXR-emitting 
energy  is  nominal  (Figure  1(a)).  Shih  et  al.  (2009)  show  that 
this  flare  has  a  ratio  of  accelerated  >20  MeV  nucleon-1 
ions  to  accelerated  relativistic  electrons  that  is  comparable 
to  other  gamma-ray  flares,  suggesting  that  this  flare  is  an 
outlier  because  it  has  a  relative  deficiency  in  higher-energy 
particles,  both  ions  and  electrons.  In  particular,  the  stated 
ion  energy  content  may  be  a  significant  underestimate  if  the 
ion  spectrum  is  steeper  than  the  S  —  4  power  law  assumed 
in  the  extrapolation  from  the  >30  MeV  proton  energy  value 
that  is  obtained  from  the  observed  neutron-capture  line 
fluence. 

5.  Event  12.  This  X17  event,  on  2003  October  28,  has 
relatively  high  CME  and  magnetic  energies  (Figure  3(a); 
Table  1).  It  was  one  of  three  events  from  the  “Halloween” 
active  region  of  2003  October-November,  in  which  a  very 
high  non-potential  magnetic  energy  value  (4  x  1033  erg; 
see  Figure  2(d))  was  inferred.  This  event  is  therefore  an 
“outlier”  simply  because  it  was  a  very  large  event;  there  are 
no  particularly  unusual  ratios  of  energetic  components. 

6.  Event  17.  The  data  set  for  this  XI. 6  event,  on  2004  July  15, 
is  quite  extensive,  with  all  components  measured  except 
CME  and  SEP  energies  (which  are  understandably  absent 
given  the  ~E45°  location  of  the  flare).  The  event  has  a 
slightly  low  ratio  of  bolometric  radiant  energy  to  magnetic 
energy  (Figure  2(d)),  and  slightly  high  ratios  of  flare- 
accelerated  particle  energy  to  bolometric  radiated  energy 
(Figure  3(c))  and  energy  radiated  by  SXR-emitting  plasma 
to  bolometric  radiated  energy  (Figure  3(d)).  It  has  the  equal- 
lowest  bolometrically  radiated  energy  content,  and  the 
lowest  thermal  energy  content  of  SXR-emitting  plasma,  of 
any  event  studied  (Table  1),  but  interestingly  does  not  show 
as  an  outlier  in  any  other  plots  in  which  it  appears,  notably 
the  plot  of  thermal  energy  content  versus  thermal  radiated 
energy  (Figure  1(b))  and  the  plot  of  thermal  energy  content 
versus  accelerated  particle  energy  (Figure  1(c))  (although 
it  is  barely  inside  the  2a  ellipse  in  both  of  these  plots).  As 
with  Event  8,  these  event  properties  may  be  explained  by  the 
relatively  short  SXR  duration  of  this  flare — only  4  minutes 
from  the  GOES  1-8  A  peak  to  50%  of  peak  flux.  Although 
the  peak  plasma  temperature  derived  from  the  GOES  fluxes 
was  22  MK,  the  temperature  stayed  above  1 5  MK  for  only 
about  8  minutes.  The  short  duration  and  low  temperatures 
lead  the  F1SM  modeling  process  (Section  2.2)  to  assign 
a  correspondingly  low  estimate  of  the  bolometric  radiant 
energy  (Chamberlin  et  al.  2012). 

7.  Event  18.  Data  for  this  M7.1  event,  on  2004  July  25,  is 
rather  limited  (Table  1).  The  event  appears  as  an  outlier 
in  Figure  2(d),  due  to  a  low  ratio  of  bolometric  energy  to 
magnetic  energy.  Table  1  shows  that  this  event  had  one  of 
the  highest  inferred  non-potential  magnetic  energy  contents 
of  all  the  events  studied,  while  Figure  1(a)  shows  very  small 
values  of  both  radiation  in  the  GOES  1-8  A  waveband  and 
total  energy  radiated  by  the  SXR-emitting  plasma  (although 
the  ratio  of  the  two  is  nominal).  RHESSI  did  show  a  flare  flag 
starting  at  05:37  UT,  but  entered  night  at  05:42  UT,  three 


minutes  after  the  GOES  start  time;  therefore  no  reliable 
electron  or  ion  energy  measurements  are  available.  Further, 
no  CME  was  observed  for  this  event,  which  was  located  at 
approximately  W30°  longitude.  As  a  result,  it  is  not  known 
whether  the  event  was  simply  very  weak  radiatively,  or 
whether  only  a  relatively  small  part  of  the  non-potential 
magnetic  energy  available  was  released. 

8.  Event  25.  For  this  X7.1  event,  on  2005  January  20,  the 
ratios  of  nonthermal  particle  energy  to  bolometric  radiated 
energy  (Figure  3(c)),  SXR-radiated  to  peak  thermal  energy 
(Figure  1(b)),  and  SXR-radiated  to  nonthermal  particles 
(Figure  1(d))  are  nominal.  However,  the  ion-to-electron 
energy  ratio  (~5:1;  Figure  2(a);  Table  1)  is  conspicuously 
high;  indeed,  this  event  had  one  of  the  highest  ion  energy 
contents  measured.  This  high  ion-to-electron  energy  ratio 
may  in  part  be  due  to  the  presence  of  high-energy  protons 
(inferred  from  the  2.223  MeV  line)  for  several  minutes  after 
most  of  the  electron-ion  bremsstrahlung  had  dissipated. 
The  event  also  has  a  very  low  ratio  of  CME  energy  to 
magnetic  energy  (Figure  3(a)),  but  a  very  high  ratio  of 
SEP  energy  to  CME  energy  (Figure  2(b)).  The  ratio  of 
bolometric  radiated  energy  to  free  magnetic  energy  is 
nominal  (Figure  2(d)),  but  it  is  the  only  event  with  a  ratio  of 
bolometric  radiated  energy  to  CME  energy  that  is  greater 
than  100%  (Figure  3(b)).  Together,  these  results  show  that 
the  main  component  that  makes  this  event  an  “outlier”  is 
the  low  CME  energy.  Note,  however,  that  the  CME  kinetic 
energy  in  this  event  is  very  uncertain. 

The  reasons  for  these  unusual  circumstances  are  not  com¬ 
pletely  understood.  Much  has  been  written  about  this  event 
(e.g.,  Grechnev  et  al.  2008,  and  references  therein)  but 
the  acceleration  of  the  intense  flux  of  SEPs  has  still  not 
been  fully  resolved.  It  has  been  suggested  that  in  this 
event  the  SEPs  were  accelerated  at  the  flare  site  rather  than 
in  the  CME  shock  (Lin  2005a,  2005b).  In  similar  vein, 
Simnett  (2006,  2007)  concluded  that,  “the  relativistic  pro¬ 
tons  were  not  accelerated  by  the  CME-driven  shock.”  The 
event  produced  a  cosmic-ray  ground-level  enhancement 
that  is  among  the  largest  recorded  in  the  history  of  cosmic- 
ray  measurements  (Mishev  et  al.  2011).  However,  based  on 
the  particle  spectrum  measured  over  a  wide  rigidity  range 
(1-20  GV),  Morgan  &  Lopate  (2008)  state  that,  “The  2005 
January  20  GLE  was  an  unusual  event  in  its  intensity  and 
brevity,  placing  it  on  the  outer  edges  of  parameter  space 
for  shock  acceleration  to  GeV  energies,  but  still  not  requir¬ 
ing  a  different  process,  i.e.,  direct  solar-flare  acceleration.” 
Similarly,  Reames  (2009)  finds  that  SEPs  in  this  event,  as 
in  other  ground-level  events  (GLEs),  are  produced  after  the 
shock  onset  and  do  not  require  a  separate  non-shock  in¬ 
jection.  On  the  other  hand,  Moraal  &  McCracken  (2012) 
conclude  that  neutron  monitor  observations  of  this  event 
indicate  two  separate  pulses  of  high-energy  particles,  one 
accelerated  by  the  flare  and  a  second  accelerated  by  the 
CME-driven  shock. 

9.  Event  30.  This  M6.4  event,  on  2005  August  25,  had  a 
very  low  peak  thermal  energy  content  and  a  relatively  low 
ratio  of  the  total  SXR-emitted  energy  to  the  peak  thermal 
energy  content  (Figure  1(b)).  It  also  has  a  low  ratio  of 
total  SXR-emitted  energy  to  nonthermal  particle  energy 
(Figure  1(d)).  No  CME  or  SEP  data  are  available  for  this 
east-limb  event.  For  this  event,  the  GOES  flux  decayed  from 
the  M6.4  peak  to  the  C6  level  (10%  of  the  peak  flux)  in  only 
~  1 4  minutes,  so  that  the  low  values  of  the  radiated  energy  in 
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Figure  4.  Bar  chart  showing  the  (logarithmic)  average  energies  of  the  different  components  for  the  six  events  for  which  values  were  obtained  for  all  components — Events 
13,  14,  20,  23,  25,  and  38.  The  short  thin  bars  show  the  ±  la  logarithmic  scatter  of  the  energies  of  the  six  events. 


SXR-emitting  plasma  are  due  to  this  simply  being  a 
relatively  short-lived  event. 

6.  CONCLUSIONS 

Despite  the  rather  large  uncertainties  in  the  individual  mea¬ 
surements  used  in  this  analysis,  the  relatively  large  number  (38) 
of  events  nevertheless  allows  us  to  reach  some  general  conclu¬ 
sions  about  the  “typical”  ratios  of  various  energetic  components 
in  large  SEEs.  We  have  found  the  following  general  statements 
to  hold. 

1.  Figure  1(b)  shows  that  the  total  energy  radiated  by  the 
SXR-emitting  plasma  over  the  course  of  the  event  exceeds, 
by  about  half  an  order  of  magnitude  (R  —  2.8  in  Table  3), 
the  peak  energy  content  of  the  thermal  plasma  that  produces 
this  radiation.  This  reinforces  the  conclusions  of  Moore 
et  al.  (1980)  that  some  form  of  energy  is  continuously 
supplied  to  this  hot  plasma  throughout  the  event. 

2.  Figures  1(d)  and  3(c)  show  that  the  energy  content  in 
flare-accelerated  particles  (electrons  and  ions)  is  suffi¬ 
cient  to  create  not  only  the  total  energy  radiated  by  the 
SXR-emitting  plasma,  but  also  the  total  bolometric  radi¬ 
ated  energy  of  the  event. 

3.  Figure  2(a)  shows  that  the  energy  contents  of  flare- 
accelerated  ions  and  electrons  are  comparable  at  the  order- 
of-magnitude  level.  This  result  supports  the  earlier  claims  of 
Ramaty  et  al.  (1995)  and  Ramaty  &  Mandzhavidze  (2000) 
and  has  significant  consequences  for  acceleration  models. 

4.  Figure  2(b)  shows  that  the  SEP  energy  is  typically  a  few 
percent  (R  —  0.04  in  Table  3)  of  the  CME  kinetic  energy 
in  the  solar  wind  rest  frame,  a  result  with  implications  for 
shock-acceleration  models  of  interplanetary  particles. 

5.  The  combination  of  Figures  2(d),  3(a),  and  3(c)  shows 
that  the  available  magnetic  energy  is  indeed  sufficient 
to  power  the  thermal  plasma,  flare-accelerated  particles, 
and  the  CME.  Although  some  “double-counting”  may  be 
involved  in  summing  these  energy  components  (e.g.,  both 
the  flare-accelerated  particles  and  the  CME  may  transfer 
energy  to  the  ambient  plasma;  see  Emslie  et  al.  2005),  this 
result  nevertheless  conforms  to  the  widely  held  view  that 


the  source  of  the  energy  released  in  SEEs  lies  in  stressed 
magnetic  fields. 

Figure  4  shows  the  logarithmic  average  (i.e.,  the  geometric 
mean)  of  the  energies  with  ±ler  logarithmic  scatter  of  the 
various  energy  components  for  the  six  events  (Events  13,  14, 
20,  23,  25,  and  38)  for  which  all  energetic  components  were 
measured.  (Events  6  and  12  were  not  included,  since  for  these 
events  some  of  the  components  were  determined  only  as  upper 
or  lower  limits — see  Table  1.)  This  figure,  coupled  with  the 
overall  ratios  summarized  in  Table  3,  succinctly  demonstrates 
how,  in  very  approximate  terms,  the  available  magnetic  energy 
gets  distributed  in  a  “typical”  flare  in  our  sample: 

1.  Of  the  ~1033  erg  of  available  non-potential  magnetic 
energy,  approximately  30%  is  released  in  the  SEE,  with  the 
remainder  staying  in  the  active  region  as  stored  magnetic 
energy.  Of  the  ~30%  that  is  released,  some  80%  (~25% 
of  the  available  energy)  is  released  in  the  CME  (mostly 
as  kinetic  energy)  and  approximately  20%  (~5%  of  the 
available  energy)  is  released  as  flare-accelerated  particles, 
roughly  evenly  distributed  between  electrons  and  ions. 

2.  All  of  the  energy  in  the  flare-accelerated  particles  appears  to 
ultimately  emerge  as  radiation  across  a  wide  range  of  wave¬ 
lengths,  from  optical  to  SXRs  (Emslie  et  al.  2005).  How¬ 
ever,  only  about  one-third  of  the  energy  in  flare-accelerated 
particles  (~2%  of  the  available  stored  energy)  is  ultimately 
radiated  from  high-temperature  SXR-emitting  plasma.  The 
maximum  amount  of  energy  stored  as  enhanced  thermal  en¬ 
ergy  in  the  SXR-emitting  plasma  is  ~  1  %  of  that  released, 
and  the  amount  of  energy  radiated  in  the  diagnostic  GOES 
1-8  A  waveband  is  only  about  5%  of  the  total  energy  radi¬ 
ated  by  the  SXR-emitting  plasma,  or  ~0. 1  %  of  the  available 
magnetic  energy. 

3.  Because  of  the  need  for  a  CME  to  “overtake”  the  solar  wind 
and  form  a  shock  front  where  SEPs  can  be  accelerated,  only 
about  two-thirds  of  the  kinetic  energy  carried  by  the  CME 
(~15%  of  the  available  non-potential  magnetic  energy)  is 
available  for  SEP  acceleration.  The  SEP  production  process 
is  in  turn  ~4%  efficient,  so  that  only  about  half  a  percent 
of  the  released  magnetic  energy  ultimately  appears  in  the 
form  of  SEPs. 
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Although  for  completeness  we  have  listed  Spearman’s  rank 
correlation  coefficient  (p)  values  in  Table  3,  little  significance 
can  be  attached  to  these  values  other  than  for  the  obviously 
tight  correlations  between  parameters  that  are  essentially  de¬ 
rived  from  the  same  data  (e.g.,  the  bolometric  emission,  the 
energy  radiated  from  the  SXR-emitting  plasma,  and  the  energy 
radiated  in  the  GOES  1-8  A  band,  all  of  which  are  dependent 
on  GOES  SXR  flux  measurements).  Any  correlations  among 
independent  components  are  masked  by  the  large  uncertainties 
in  the  individual  measurements  used  in  the  various  scatter  plots. 
Progress  in  this  direction  will  require  sampling  of  events  over  a 
much  larger  range  in  flare  size  to  determine  if  the  distribution  of 
energies  among  the  different  components  found  here  for  large 
events  is  preserved  for  smaller  events,  as  would  be  expected 
for  a  “ big  SEE  syndrome .”  Such  a  project  is  the  next  step  to¬ 
ward  a  more  comprehensive  understanding  of  energy  release 
in  SEEs. 

Using  data  from  the  new,  more  sophisticated,  instruments  that 
are  now  available  will  allow  more  accurate  energy  estimates  to 
be  made  of  some  of  the  components.  For  example,  the  detailed 
differential  emission  measure  analysis  now  possible  using  data 
from  the  EUV  Imaging  Spectrometer  (Culhane  et  al.  2007)  on 
Hinode,  and  from  the  Atmospheric  Imaging  Assembly  (Lemen 
et  al.  2012)  and  the  Extreme  Ultraviolet  Variability  Experiment 
(Woods  et  al.  2012)  on  the  SDO,  should  lead  to  better  estimates 
of  the  energies  in  the  thermal  plasma.  Vector  magnetograms 
from  the  HMI  (Scherrer  et  al.  2012)  on  SDO  allow  for  more 
accurate  estimates  of  the  energy  in  the  non-potential  magnetic 
field. 

Other  energy  components  not  considered  in  our  analysis  may 
be  found  to  contain  significant  total  energy  and  should  be  in¬ 
cluded  in  any  future  compilation  of  global  energetics.  These 
include  the  turbulent  mass  motions  revealed  by  the  broaden¬ 
ing  of  atomic  lines  seen  with  EUV  and  X-ray  spectrometers 
(e.g.,  Phillips  et  al.  2008),  and  the  cumulative  heating  of  CME 
plasma  reported  by  Murphy  et  al.  (2011)  to  be  comparable  to 
(or  even  greater  than)  the  CME  kinetic  energy.  Another  aspect 
not  discussed  here  is  the  question  of  a  second  flare  phase  that, 
according  to  Woods  et  al.  (2011)  and  Su  et  al.  (2012),  can 
release  a  similar  amount  of  energy  as  in  the  initial  phase.  Never¬ 
theless,  we  believe  that  the  order-of-magnitude  comparisons  of 
energetic  components  presented  herein  represent  a  significant 
advance  in  our  understanding  of  the  nature  of  energy  release 
in  SEEs. 
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both  perpendicular  and  parallel  to  the  lines  of  constant  ratio,  and 
the  Spearman’s  rank  correlation  coefficient  p. 

The  coordinates  of  the  logarithmic  centroid  in  each  plot  are 
given  by 


l°Sio  -^centroid 


log io  ^centroid 


1 

-£log10X;, 

i=l 

l  N 

*I>*o*. 

i=i 
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where  N  is  the  number  of  events  for  which  there  are  viable 
measures  of  both  components  included  in  the  scatter  plot  in 
question. 

The  lines  of  constant  ratios  ( R  =  100%,  10%,  and  1%) 
between  the  X  and  Y  components  satisfy  the  following  relation: 


R  =  Y/X 

or  log10  R  =  log10(T/ X) 

=  log  10  Y  -  log10  X.  (A2) 

Lines  of  constant  logarithmic  average  event  energy  (log10  A), 
shown  in  the  plots  as  tick  marks  every  order  of  magnitude  along 
the  lines  of  constant  ratio,  are  defined  as  follows: 


logio  A  =  (log10  X  +  log10  Y)/2 

=  log10  \[XY .  (A3) 

The  rms  deviations  of  the  points  perpendicular  to  and  parallel 
to  the  line  of  constant  ratio  passing  through  the  centroid  are 
defined  by 


rmsi  =  (1/A0  /y:  (log10  R;  -  log10  Rcentroid)2 

=  (1/A) Oogi0(F,/ X,)  -  log10(Tcentroid/Xcentroid))2 


(A4) 


and 


rms ||  —  (1/A)  / 'y  '  (logio  A;  logio  ^  centroid)^ 

=  (1/2A)  [logical  U)  -  log 1 0 (Xcentroid ^centroid ) ] 2 • 


(A5) 


These  rms  values  were  used  to  draw  the  ellipse  with  axes  of 
2  x  rmsx  and  2  x  rmsy  around  the  logarithmic  centroids  in  each 
plot. 

The  Spearman’s  rank  correlation  coefficient  p  is  calculated 
by  first  assigning  ranks  jc,-  (=  1 , . . . ,  A)  and  yt  (=  1 , . . . ,  N)  to 
the  X  and  Y  values,  respectively,  of  the  Appoints  used  in  the  plot 
in  question.  The  ranks  are  assigned  such  that  X (xj )  X(x,+| ) 
and  Y(yt)  ^  T(y;+1),  with  “ties”  assigned  the  average  rank  of 
the  tied  values.  (Note  that  the  rank  order  does  not  depend  on 
whether  we  use  the  X,  values  or  their  logarithms  log10  X,-.)  The 
Spearman’s  rank  correlation  coefficient  is  then  calculated  as  the 
correlation  coefficient  of  the  ranks: 


lEife  ~  -t)(y/  -  y) 

yELfc-*)2Ef=i(v;-y)2’ 


(A6) 


Here,  for  definiteness,  we  provide  the  equations  used  to 
determine  the  values  listed  in  Table  3:  the  (logarithmic)  centroid 
energies  XcentIOi(i  and  Tcentroid>  their  ratio  R.  the  rms  (la)  values 
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where  x  and  y  are  the  means  of  the  ranks  x,  and  y,,  respectively. 
For  a  monotonic  dependence  of  Y  on  X,  p  =  1,  even  if  the 
variables  do  not  obey  a  perfect  linear  correlation. 
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